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Abstract 

We detail the automatic construction of R matrices corresponding to (the 
tensor products of) the (0 m \6t n ) families of highest-weight representations 
of the quantum superalgebras U q [gl(m\n)]. These representations are ir- 
reducible, contain a free complex parameter a, and are 2 mn dimensional. 
Our R matrices are actually (sparse) rank 4 tensors, containing a total of 
2 4mn components, each of which is in general an algebraic expression in 
the two complex variables q and a. 

Although the constructions are straightforward, we describe them in 
full here, to fill a perceived gap in the literature. As the algorithms are 
generally impracticable for manual calculation, we have implemented the 
entire process in Mathematica; illustrating our results with (7 9 [(7Z(3|1)]. 



1 Introduction 

Broadly, R matrices are solutions to the various versions of the Yang-Baxter 
equation, and as such, are of great interest in mathematical physics and knot 
theory (see, e.g. [@), both in their algebraic (i.e. "universal") forms, and 
in their (matrix) representations (i.e. "quantum" forms), useful for explicit 
computations. Here, we will be specifically concerned with quantum R matrices 
associated with the quantum superalgebras U q [gl(m\n)]. 

Although much is known about the origin and properties of quantum super- 
algebra R matrices (e.g. |l9[ ] provides universal R matrices), explicit examples of 
their quantum R matrices are rare in the literature, due largely to the computa- 
tional effort involved in obtaining them. This paper describes the automation of 
an algorithm to generate a suite of explicit quantum R matrices for U q [gl(m\n)]. 
As readers of this organ may not be familiar with these algebraic structures, we 
provide a full description of their details. 
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Specifically, we construct trigonometric R matrices R m ' n (u) corresponding 
to the a-parametric highest weight minimal representations labeled (0 m |d„) of 
the U q [gl(m\n)]. These irreducible representations are 2 mn dimensional, and 
contain free complex parameters q and a; the real variable u is a 'spectral' pa- 
rameter. Quantum R matrices R m ' n are immediately obtainable as the spectral 
limits u — ► oo of the R m,n {u). 

Our R matrices are in fact graded, as they are based on graded vector spaces, 
hence they actually satisfy graded Yang-Baxter equations. However, it is a 
simple matter to remove this grading and transform them into objects that 
satisfy the usual Yang-Baxter equations. 

The constructions have been implemented in Mathematica, and results ob- 
tained for n = 1 and m = 1,2, 3, 4; we illustrate the algorithms using C/ q [gZ(3|l)]. 
Full listings of all our R matrices have been announced in . 

As they are solutions to Yang-Baxter equations, our R matrices are of im- 
mediate practical interest. Firstly, they are of physical interest in that they are 
applicable to the construction of exactly solvable models of interacting fermions. 
Corresponding to R m,1 (u), we may construct an integrable 2 m state fermionic 
model on a lattice. Models associated with J7 g [gZ(2|l)] and J7 g [pi(3|l)] have been 
discussed in Q and respectively. The £/ g [gZ(4|l)] case has an elegant inter- 
pretation in terms of a 2-leg ladder model for interacting electrons: a discussion 
of this is provided in [|| . 

Furthermore, corresponding to each R m,n , we may obtain a polynomial 
'Links-Gould' link invariant LG m,n plL cf. the celebrated Jones polynomial. 
These LG mn are two- variable, integer-coefficient Laurent polynomials, and are 
generally substantially more powerful than the Jones polynomial in distinguish- 
ing knots. (LG 1,1 degenerates to the well-known Alexander-Conway polynomial 
in the single variable q 2a (cf. 0).) A fuller documentation of the suite LG m,n 
has been provided by myself in collaboration with Louis Kauffman and Jon Links 
in ||, |7], [|. Although the LG m ' n are far from being complete invariants, as 
they can distinguish neither mutants nor inversion [|, it turns out that even 
LG 2:1 is in fact more powerful than the well-known two-variable HOMFLY and 
Kauffman invariants, being able to distinguish (including chirality) all prime 
knots of up to 10 crossings ||. Their evaluation also involves automatic sym- 
bolic computation, but the computational aspects are comparatively pedestrian. 

Lastly, we mention explicitly that this paper contains no new theorems, 
although it does contain two new technical lemmas, proven in Appendix |a|. It 
is primarily intended to provide a proper foundation for the results presented 
in [^,0, although it also serves as a tutorial on an application of symbolic 
computation. Whilst it specifically pertains to representations of U q [gl(m\n)], 
many of the algorithms have a much broader application. 

The following subsections provide a synopsis of the paper. 
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1.1 Algebraic overview 

Fixing to and 77, we are initially interested in a 2 mn dimensional vector space 

V that is a module for the U q [gl(m\n)] minimal typical highest weight repre- 
sentation A = (0 m |d„). The algebra contains a free complex variable q, whilst 
the representation tt\ acting on V contains a free complex variable a. Our 

V is actually (Z2) graded; this ensures compatibility with the (Z2) grading of 
U q [gl(m\n)]. 

Using the properties of U q \gl(m\n)], we apply a version of the Kac induced 
module construction (KIMC) [l7| to establish a (weight) basis {|£)}f=i f° r 
V. This involves postulating |1) as a highest weight vector, and recursively 
acting on |1) with all possible distinct products of simple lowering generators 
E a+1 a to define the other basis vectors, normalising as we go. This construction 
requires a 'Poincare-Birkhoff-Witt (PBW) lemma' for U q [gl{m\n)\ || |28 



i.e. 



a set of commutations sufficient to transform any product of algebra generators 



into a normal form (see §4.2), together with a statement that the algebra is 
spanned by the set of all such normal forms. 

Where V has a graded weight basis {|«)}| =1 , the tensor product module 
V ® V has a natural 2 2mn dimensional basis {\i) ® |j)}f J=1 , which inherits 
a weight system and a grading from V. For our particular representation, the 
orthogonal decomposition of V ®V is known JllJ, and contains no multiplicities, 
viz: 

v®v = ©v fc , 

fe 

where the submodule Vk has highest weight A^, and these Xk are known, and 
all distinct. To build R matrices acting on V (g> V, we require an alternative, 
orthonormal weight basis 03 = {J k *Bfc for V (8) V, corresponding to this decom- 
position, viz Q5fe is a basis for Vk- Again using the KIMC, the basis vectors 
of each 55fc are derived as linear combinations of the form 0ij(\i) Cg) \ j)), where 
the coefficients Oij are algebraic expressions in q and a. This process initially 
yields a basis 25fc that is not necessarily orthonormal, so we also apply a Gram- 
Schmidt process to orthonormalise 25^ into Q3fe. The desired R matrix is then 
a weighted sum of projectors onto these Vk, where the weights are eigenvalues 
of the appropriate second order Casimir invariants. 

The algebraic structure of U q [gl(m\n)} is detailed in §pj, and an introduction 
to its highest weight representations is provided in §^|. In we provide a normal 
ordering and a PBW lemma for U q [gl(m\n)]. The construction of our particular 
(6 m |d„) representations is detailed in §|B|. §|| describes the construction of the 
bases 03^, and describes the construction of projectors and R matrices. 



3 



1.2 Implementation and results 

Explicit computations within the representation theory of quantum superalge- 
bras are tedious and error-prone when performed manually. The dimensions 
of representations are generally large, and in our case, we have the presence of 
the two variables q and a; these generally manifest themselves in complicated 
rational algebraic expressions, whose symmetries must be continually identified 
and exploited to avoid the arising of intractable messes of algebra. 

The construction of the basis {|i)}f =1 involves many applications of the 
PBW lemma to simplify long strings of algebra generators. This is compu- 
tationally expensive; firstly as the simplification involves a minimally-efficient 
sorting process, and second as it involves a geometric explosion in the number 
of terms being sorted. 

The construction of the weight space bases 5^ is nontrivial, as each basis 
vector of each <Bfc generally contains many terms of the form Oij(\i) <E> | j))> where 
the coefficients % are generally complicated rational algebraic expressions in q 
and a. (That said, we have avoided the more theoretically difficult situation of 
computing weight space decompositions in cases where there are weight multi- 
plicities in the underlying carrier space V.) Although the R matrices have 2 4m ™ 
components, Nature is kind to us in that most of these components are zero, 
and those that are not are generally simpler than the 6%j . 

To the best of our knowledge, computer implementation of the algebraic 
structures and algorithms described herein has not previously been achieved. 
We have implemented the entire process as a suite of Mathematica functions; 
the thousands of lines of code perform algebraic computations that a human 
being could not ever realistically expect to perform correctly. 



From §2^3 onwards, we use [/ g [gZ(3|l)] to illustrate our results. These are 
summarised in Appendix [b|; where we list the explicit matrix elements for the 
generators of the underlying 8 dimensional representation, orthonormal bases 

for the 4 submodules V& C V <S> V, the components of the associated 4 pro- 
jectors Pk onto the and finally, the trigonometric and quantum R matrices, 
R 3:1 (u) and R 3,1 , respectively. 

Whilst there are no theoretical limits to m and n, a current practical limit for 
computation is ran ^ 4. This is convenient, as an immediate application || of 
the material critically requires J? 4,1 (u). Although translation of the interpreted 
Mathematica code into a compiled language would increase the speed of the 
computations, storage requirements would still limit ran to perhaps 7 in the 
general case. 

Further discussion of implemcntational issues and results is provided in §pl 
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2 The quantum superalgebras U q [gl(m\n 



The algebraic structures labeled U q [gl(m\n)] are quantum superalgebras^ de- 
scribed in many places, e.g. |24], |2f| |2(|, and in the book (3[ see §6.5]. 

For our purposes, to and n are positive integers, to be regarded as fixed, 
and q is to be regarded as a nonzero compl ex v ariable. As U q [gl{m\n)\ may be 
unf ami liar to the readers of this organ, in §2J we introduce its phylogeny, and 
in §^2, we provide a full d escr iption of its structure in terms of ge nera tors and 
relations. Beyond that, in §2^3 we describe its root system, and in §2^ we show 
how it may be regarded as a Hopf (super) algebra. 



2.1 The phylogeny of U q [gl{m\n)] 

1. Where n is a positive integer, recall that the Lie algebra gl{n) is equivalent 
to the usual (complex) vector space of n x n (complex) matrices augmented 
by a 'vector multiplication' operation which is the usual matrix multipli- 
cation. gl{n) is of course a unital algebra, and is of dimension n 2 and 
rank n — 1. The n 2 generators {e a t,}™ b=1 of gl(n) satisfy a commutation 
relation: 

[e\,e c d }=S c b e a d -S a d e c b , 

where [•, •] is the usual commutator (bracket), defined for X, Y 6 gl{n) by: 

[X, Y] = XY - YX. 



2. Letting both to and n be positive integers, the Lie superalgebra gl(m\n) 
may be obtained from gl(m + n) by retaining the generators {e a &}™ ; ^™ 1 , 
but modifying the definition of the commutator bracket and commutation 
relations to include some 'parity factors' of ±1. Specifically, we have the 
commutation relation: 

[e\,e c d ]=8 c b e a d -(-) [eab][eCd] S a d e c b , (1) 



where [•, •] is now the graded commutator (bracket), defined for homoge- 
neous (see below) X,Y £ gl(m\n) by: 

[X,Y] = XY- (-)Ml y 'FX, (2) 

and extended by linearity. In both (|l|) and (||), [X] e {0, 1} refers to the 
grading of the homogeneous element X . For this reason, Lie superalgebras 
are sometimes called "graded Lie algebras". From the gl(m + n) case, we 
see that gl(m\n) is of dimension (to + n) 2 and rank m + n — 1. 



1 These structures are sometimes called "quantum supergroups", but they are actually 
(associative, noncommutative) algebras. 
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3. U[gl(m\n)] is then the usual universal enveloping algebra obtained from 
gl(m\n) by regarding the gl(m\n) generators as letters in words contained 
in U[gl(m\n)], where the (graded) commutator bracket acts as a relation to 
reduce the algebra somewhat. U[gl(m\n)] is infinite dimensional, although 
of finite rank, viz, again m + n — 1. 

4. The quantum superalgebra U q [gl(m\n)] is then a so-called 'g-deformation'Q 
of U[gl(rn\n)], which maintains its viability as a Hopf (super)algebra struc- 
ture (see below) JT^ . Roughly speaking, the deformation amounts to 
'exponentiation by q'; indeed U[gl(m\n)] may be recovered as the limit 
q — > 1 of U q [gl(m\n)]. U q [gl(m\n)] is of course also infinite dimensional, 
and again of rank m + n — 1 . 



2.2 Generators and relations for U q [gl(m\n)] 

Following Zhang (2^, ppl237-1238], we provide a full description of U q [gl(m\n)] 
in terms of generators and relations. For various invertible X , we will repeatedly 
use the notation X = X~ l . 

2.2.1 U q [gl(m\n)] generators 

Where X = {1, ...,m + n} is the set of the gl(m\n) indices, we define a Z 2 
grading [•] : I -> Z 2 : 

, , A J if a ^ m even indices 

[1 else odd indices. 

Throughout, we shall use dummy indices a, b G X where meaningful. A set of 
(m + n) generators for U q [gl(m\n)] is then: 

K a , m + n 

E b a , a<b i(m + n)(m + n 

E a b, a<b i(m + n)(m + n 

Let us now introduce the notation, for any a G I 



For any power AT, replacing q with q w immediately shows that (q a ) = (q )oj 
so we may write with impunity, specifically, we will write <f a = g^ 1 . Next, 
an equivalent notation for K a is q® a ; where the exponential is defined in the 
usual manner as an infinite sum, thus powers are meaningful; specifically, 
we will often be working with N G |Z. Thus, under the mapping q >-^q, K a is 
mapped to K a , where we intend K a = K~ x . As expected, for arbitrary powers 
M, N, we have: 

K?K? = K?+« where K° a = Id, 

where Id is the C/ 9 [c//(m|n)] identity element. Apart from N G N, powers (i.e. 
products) of the non-Cartan generators (E a i,) N for a ^ b, are not meaningful. 

2 We might say 'quantum deformation' heregbut the relation to quantum mechanics is more 
of analogy than of rigor. 



Cartan 
1) lowering ^ . (3) 
1) raising 



On the generators we define a natural Z2 grading in terms of the grading on 
the indices: 



K] = 0, [E\}^[a} + [b] (mod 2), (4) 

where the former may be seen as a special case of the latter by setting a = b and 
making the identification K a = g(~) w - B °< i . We use the terms "even" and "odd" 
for generators in the same manner as we do for indices. Elements of U q [gl(m\n)] 
are said to be homogeneous if they are linear combinations of generators of 
the same grading. The product XY of homogeneous X, Y £ U q [gl(m\n)] has 
grading: 

[XY]±[X] + [Y] (mod 2). (5) 

Thus, for example, inspection of (|4|) and (|J) shows that we may cheerfully 
substitute [E a c ] for [E a b E b c ]. Further, we also have the following useful results 
for a < b < c < d: 

[E a d ] [E\] = [E\] and [E\] [E\] = [E\] [E c d ] = 0. 



2.2.2 U q [gl(m\n)] simple generators 

The full set of generators (||) includes some redundancy; some can be regarded 
as simple in that the rest may be expressed in terms of them. We shall call the 
following subset of 3(to + n) — 2 generators the U q [gl(m\n)] simple generators: 

K a , m + n Cartan ) 

E a+1 a , a < m + n m + n — 1 simple lowering > . 

E a a+ i, a < m + n m + n — 1 simple raising J 

The fact that there are m + n — 1 simple lowering generators indicates that 
U q [gl(m\n)] has rank m + n — 1. Note that there are only two odd simple 
generators: E m+1 m (lowering) and E m m+1 (raising). 



2.2.3 U q [gl(m\n)] nonsimple generators 

In the gl(m\n) case, the remaining nonsimple (non-Cartan) generators satisfy 
the same commutation relations as the simple generators. The situation is 
very different for U q [gl(m\n)]; the nonsimple generators do not satisfy the same 
commutation relations as do the simple generators. Instead, they are recursively 
defined in terms of sums of products of the simple generators (see |27], pl971, (3)] 
and |^||, pl238, (2)]). Strictly speaking, they are not explicitly required for the 
definition of the algebra; their use can help simplify otherwise large expressions. 

To whit, a set of U q [gl(m\n)] nonsimple generators may be defined recursively 
for a < b by: 



(a) E h a = E b c E c a — q c E c a E b c nonsimple lowering 

(b) E a t, — E a c E c b — q c E c bE a c nonsimple raising, 




where a < c < b; viz c is an arbitrary index, we do not intend a sum here. 
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In we will have use for an alternative set of nonsimple generators, again 
defined recursively for a < b by: 

(a) E h Q = E 6 c E c a — q c E c a ~E b c alternative nonsimple lowering 

(b) E a 6 = E a c E c b — g c E c f,E a c alternative nonsimple raising, 

where we intend E a b = E a b when E a b is simple, viz, for any a — b\ = 1. Note 
that we use a boldface E a b where the original source |^7j uses an overline if;,; 
the use of the boldface notation saves the overline for indicating inverses. 

These definitions may be written more concisely with some more notation. 
Writing S% = sign(a — b), we may replace and (0), for all a ^ b, by: 

(a) E\ = E a c E c b -q?E c b E a c 

(b) E% = E a c E c h -gf E c 6 E a c . 

The two different sets of generators are in fact Hermitian conjugates. For 
all meaningful indices a, 6, we have: 

(£%)t = E b a , (E a 6 )t = E» ai (K^ = K», 

and these definitions ensure that {X<y = X for all U q [gl(m\n)] generators X. 
Note that these are ordinary, not Z2 graded Hermitian conjugates, meaning that 
we have (XY)i = Y^X^; expressly not (17)+ = (-)™ Y ^X^. 

Lastly, we mention a result of Zhang [P7[ Lemma 3], which gives us a more 
efficient formula than (|^b) for expanding the alternative nonsimple generators: 

E a fc = E\ + 5 h a A c ^ c h E a c , (9) 

c 

for any indices a^6, where the sum is over all c strictly between a and b. (If 
I a — fo| = 1, then the sum is ignored, and the result is trivial.) 

Note that in (^|), we have introduced the following handy notation: 

A 4 q-q, A a ±q a -q a = (-) [a] ( q -q) = (-) [a] A 
A 4 (A)- 1 A a = (A a ) _1 . 

2.2.4 The graded commutator 

The graded commutator [•,•] : £/ ? [gZ(m|n)] x Uq[gl(m\n)] —> J7 g [gZ(m|7i)], is de- 
fined for homogeneous X, Y <E [/ q [gZ(m|n)] by (0), viz: 

[X,Y] ^ XY - {-) lx]lY] YX, 

and extended by linearity. For completeness, we mention that for associative su- 
peralgebras, of which U q [gl(m\n)] is certainly an example, we have the following 
useful graded commutator identities: 

(a) [XY,Z] = X[Y,Z] + {-) [Y][Z] [X,Z]Y 

(b) [X,YZ] = [X,Y]Z +{-) [x][Y] Y[X,Z\. 
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2.2.5 U q [gl(m\n)) relations 

With this notation, we have the following U q [gl(m\n)] relations: 

1. The Cartan generators all commute; for any powers M,N: 

KfK* = K*K™. (12) 

2. The Cartan generators commute with the simple raising and lowering 
generators in the following manner: 

K a E\ ±l K a = q { f £ ^ ±l] E b b±1 . (13) 



From (p. 3D, we have the following useful interchange: 

K a E\ ±l = qf b ^ ab±l) E\ ±l K a . (14) 

In Lemma || (proved in Appendix |A.l| ) , we show that ( [l4] ) may be much 
strengthened to: 

K»E\ = <£W-^E\K», (15) 

for any meaningful indices b, c (viz b < c, b > c, and even b = c), and any 
power N. 

3. The non-Cartan simple generators satisfy the following commutation re- 
lations (this is the really interesting part!): 

[E\ +1 ,E b +\] = § a K aK a+ i-K a K a+1 

Qa - q a 

Alternatively, again employing the notation of (|l^), we may write this: 
[E\ +1 , E b+ \] = 5Z(-) la] \(-) la] E a a - (-) [a+11 £ a+1 a+1 l , (17) 

L J q 

where we have introduced the g-bracket, defined for various invertible 
X 6 Uq[gl{m\n)\, including scalars (well, scalar multiples of Id): 

4 q ~ q observe that lim [XL = X. (18) 

q-q 8->i 9 

Note that in @, Zhang || replaces (-) [a] E a a - {-) [a+1] E a+1 a+1 with 
the more convenient expression h a . 
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In fact, @ generalises to a more useful result (proven in J§] ) : 



[E a b E b } = KaKb Z KaKb 
la - q a 



A a (K a K b - K a K b ), 



(19) 



viz: 



[E a b ,E\] = (-) la| (-) w E a a - (-) [f V- 



We also have, for \a — b\ > 1, the commutations: 

E a+1 a E b+1 b = E b+1 bE a+1 a and E a a+ iE b b + i = E b b+iE a a+ \. 



(20) 



4. The squares of the odd simple generators are zero: 

{E m m+l f = (E m+1 m f = 0. 

In fact, we may show that this implies that the squares of nonsimple odd 
generators are also zero: 



{E a b f = o, 



[a] ± [b]. 



(21) 



5. Lastly, we have the U q [gl(rn\n)] Serre relations; their inclusion ensures 
that the algebra is reduced enough to be simple. For a ^ nr. 

{E a+1 a ) 2 E a±1+1 a±1 - VE a+1 a E a±1 +i ±1 E a+1 a + E a±1+1 a±1 {E a+1 a f = 

ytj a + l) Hi a±l + l ~ Vii o+l* ail + l-E' a+1 + Hi a±l + l[Hj a +l) — U. 



where to save space, we have introduced the notation: V = q + q. These 
may be more succinctly expressed using nonsimple generators. Noting 
that for a ^ m, we have q a — q a +i, the above become, again for a ^ m: 



(a) E a+1 a E a+2 a 

(b) E a a J r \E a a +2 

(c) E a+1 a ^E a+1 a 

(d) E a ~ 1 a+\E a a+ i 



qa,E a a+2E a a +l 
„ rpa+1 rpa+1 

~ Co Tpa— 1 

a+l-E' a+l- 



> (22) 



10 



Alternatively, if we define a q graded commutator: 

[X, Y] q ^XY - {-)™ Y ^qYX, 

and, equivalently, a q graded commutator by replacing q with q, then ( p2| ) 
may be more elegantly expressed: 

[x,[x,y],v = o, (23) 

where the pair (X,Y) represents the four pairs (E a a+ i, E a±1 a ±i+ 1 ) and 
(E a+1 a , E a±1+1 a ±i). Equivalently, we may exchange q and q in (^3]). Note 
that in these cases, the parity factors (— )[ x l[ y l in ( |23| ) are always +1 as 
a / m in E",, + i and E a+l a , viz the graded commutators degenerate to 
ordinary commutators. 

The a ^ m Serre relations are complemented by a pair dealing with the 

case a = m: 

T^m+l LTm+2 rnm+2 

mt-J m — 1 — m—l-*- J m 

rpm rpm—l rpm — 1 zp m 

Hi m+l-C- rn+2 — ~ & m+2-& m+1, 

more succinctly expressed as: 

L-ti- mi-fr m-lj — [ m+li-^ m+2j — U, 

where, as the generators are all odd, the graded commutators are read as 
anticommutators . 

Observe that if either m or n is 1, there are actually no U q [gl(m\n)] Serre 
relations; making life a little simpler. 

These relations tell us, in principle, how to reexpress products of simple 
generators. In general, to reexpress a product containing nonsimple generators, 
those nonsimple generators must first be recursively expanded using (||a), with 
any graded commutators expanded by linearity, before the above relations can 
be invoked. 

The above description of the U q [gl(m\n)] relations should convince the gentle 
reader that U q [gl(m\n)] has a formidable structure. To facilitate examination 
of its representation theory, in §U, we will rewrite the U q [gl(m\n)} relations into 
a PBW basis formulation, which is suitable for implementation on a computer. 
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2.3 U q [gl(m\n)] root system 

We next introduce the U q [gl(m\n)] root system, which is identical to that of 
gl(m\n). We will have use for it in §3.2, §4.1 and §5.1. 



Where the gl(m\n) Cartan subalgebra is denoted by H, its dual, the gl(m\n) 
weight space H* , has a basis given by the gl(m\n) fundamental weights {e a }a£T, 
which are lists of zeros of length m + n, with a 1 in position a. The e a inherit 
a grading from that on the indices. As H and H* are dual, where {e b b } b ^x are 



the gl(m\n) Cartan generators, we have the form: e b {e a ) 
have an invariant symmetric bilinear form (•, •) : H* x H* 

(e a ,s b )^(-)^S ab , 



S ab . On H*, we 
C , defined by: 



and extended by linearity. 

Next, to each non-Cartan gl(m\n) generator e a b , there corresponds a gl(m\n 
root a ab = e a — e b , which is the weight of e a b in the adjoint representation 
For our purposes, it is convenient to bastardise the notation. Also permitting 
a = b, we will refer to a ab as the 'weight' of e a b : 



wt(e a f,) = e a - e b = a ab , 



(24) 



indicating that within a gl(m\n) weight module (see the action of a gener- 
ator e a b sends a vector of weight 7 to another of weight 7 + a ab . 

The roots inherit a grading from the indices: [a ab ] = [e a b ]. Further, we 
assign signs to them in accordance with those of these generators, viz that if 
e a b is a lowering generator, the corresponding root a ab is said to be negative, 
written a ab -< 0, and if e a b is a raising generator, then a ab is said to be positive, 
written a ab y 0. To illustrate, weights for the J7 g [gZ(3|l)] lowering generators 
are presented in Table |]. 

X wt(X) 
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0, 


0, 


-1 1 


+1) 


E\ 


( 


0. 


-1, 





+1) 


E\ 




-1, 


0, 





+1) 


E\ 


( 


0, 


-1, 


+1 1 


0) 


E\ 


(" 


-1, 


0. 




0) 


E\ 


(" 


-1, 


+1, 





0) 



odd 



even 



Table 1: Weights (all negative) of the L/q [<?Z(3|1)] lowering generators. 



Using this notation, gl(m\n) has the following simple, positive roots: 
a a ,a+i =e a — £ a+i, a = l,...,m + n — l. 

Apart from the single odd simple root a m . m +i, the simple positive roots are all 
even. (Of various choices for Lie superalgebra root systems, this distinguished 
root system is unique in containing only one odd simple root.) 

3 We apologise for the overloading of a. The notation in this subsection will go no further. 
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Then, we define Af to be the set of gl(m\n) positive roots of grading i, and 
A + to be the union of the Af, vizQ 

At = {7: [7] = * and 7 >-()}, A + = A+ U A+ 

where [7] denotes the grading of the root 7. 

In terms of these, we define the half sums of all positive even and odd gl(m\n) 
roots, and their graded sum p: 

ft = lE^ p = \ E H [7l 7= E (-)>*> 

7 eA+ 76A+ i=o,i 

viz p = po — p\. Specifically, for gl(m\ri) (and hence for U q [gl{m\nj\), we find 
@p6207]:| 

m m+n 

Po = I / ] (m — 2a + 1) g tt + ^ (m + n - 2a + 1) £ a 

a— 1 a— m-f-1 

tn m-\-n 



pi = | E n£a - 3 E m£a 

a— 1 a— m+1 

m m+n 

(m - n - 2a + l)e a + i (2to + n - 2a + 1) e a . 



2 

a=l a— m+1 



2.4 [/Jg^mlri)] as a Hopf superalgebra 

U q [gl(m\n)] may be regarded as a Hopf superalgebra when equipped with the 
following (compatible!) coproduct A, counit e and antipode S structures. The 
material is taken from |^8|, pl238], except that we have modified the definition 
of the coproduct and antipode so that they have increased symmetry. This 
material is included for completeness; in §^|, we will only have need for the 
coproduct. Uq[gl{m\n)\ is in fact a quasitriangular Hopf superalgebra, i.e. it 
possesses a (universal) R matrix. 

We first introduce some notation. A homomorphism or antihomomorphism 
H on U q [gl(m\n)] is described as (1*2 ) graded if it is compatible with the graded 
commutator, viz: 

H([X,Y]) = [H(X),H(Y)}, 

where the latter graded commutator may even exist on an ungraded space, e.g. 
C , where it is actually trivial. 

This means that a graded homomorphism H and a graded antihomomor- 
phism A must necessarily satisfy: 

H(XY) = H(X)H{Y) and A(XY) = {-) [x]lY] A(Y)A{X), 

of which only the latter varies from the usual, ungraded situation. 

4 We apologise for the overloading of A. In practice, this A will only appear with a positive 
superscript, so it is easily distinguishable. 

5 We make a correction to which app ears to cite an error reproduced several times 

before and after, e.g. it appears in Zhang pa] . To whit, the term "2m" in the formula for p 
repeatedly appears as "m". 
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2.4.1 Coproduct A 

The coproduct (a.k.a. comultiplication) is a Z 2 graded algebra homomorphism 
A : U q [gl{m\n)\ — > U q [gl(m\n)] ® U q [gl(m\n)], denned by: 



(a) A(E a +\) = E a+1 a ®KlKl +1 +KZKl +l ®E a+1 a 

(b) A(E\ +1 ) = E\ +1 ®KlKl +1 +KlKi +1 ®E a a+1 



(25) 



(c) A(K a ) = K a ®K a , 

and extended to all of U g [gl(m\n)] by: 

A(XY) = A(X)A(Y), ior all X,Y eU q [gl(m\n)}. (26) 

Observe that A preserves grading, viz that [A(X)] = [X] for homogeneous 
X e U q [gl{m\n)], where we have [X®Y} = [X] + [Y]. 
By substitution of q N for q in (p5|c), we discover: 

and, setting N = 0, hence A(Id) = Id ® Id, as expected. 



Before proceeding, we mention that our definition in (25) is only one of 
various possibilities; we have chosen it for its symmetry. In fact, in comparison 
with the literature, our A agrees with that of |l3|, and differs from that of Zhang 
[ p8| and that of my PhD thesis 0). We mention that given any coproduct, it is 
possible to write down another coproduct structure, the "opposite coproduct" : 
A T = T ■ A. Here, the twist map T is an operator on the tensor product 
U q [gl(m\n)] ® U q [gl(m\n)], defined for homogeneous X, Y £ U q [gl{m\n)\ by: 

T{X®Y) = {-) [X][Y] {Y ®X). 

More relevant to our purposes here, we may extend the expression for the 
coproduct for simple generators to that for nonsimple generators. Firstly, as 
in (^|), writing S% = sign(a — b), we may cheerfully rewrite (|2^) for the simple 
generators E a t, for any |o — b\ = 1: 

A(E\) = E\ <g> Kl Sl Kl Sl + Kl^Kf* ® E a b . (27) 

Using this notation, we prove in Lemma [3] in Appendix |A.2j the following more 
general statement, for any valid indices a, b: 

A{E%) = E a b ®K% St Kl St +Kl St K^ St ® E\ - 

1 Q a 1 CO 1 qa 1 O 



S h a ]TA c [Kf b Kf b E\ 2) E c b Kr b Kf b ) , (28) 



where the sum ranges over all c strictly between a and 6, and is simply ignored 
if \a — b\ ^ 1. Where a — b, the statement is also true; this is made clear when 
the equivalence K a — q(' ° Ea " is noted. 

Lastly, we apologise for even further overloading the definition of A; to be 
sure, the coproduct will only appear with parentheses enclosing its argument. 
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2.4.2 Counit e 



The counit e : U q [gl{m\n)\ — > C , is also a Z 2 graded algebra homomorphism, 
denned by: 

e{E a± \) = e{E a a±1 ) = 0, e{K a ) = 1, 

and extended to all of f7 g [gi(m|n)] by s(XY) = e(X)e(Y). Again, we have 
s(K^) — 1, and, setting N = 0, thus e(Id) = 1, as expected. We apologise 
for overloading the dehnition of s as the counit with the gl(m\n) fundamental 



weights (see §2.3). As we shall have no further use for the counit, we arc safe. 



2.4.3 Antipode S 

Lastly, the antipode S : U q [gl(m\n)] — > U q [gl(m\n)], is a Z 2 graded algebra 
antiautomoTpbisiai, defined by: 

S(E a+1 a ) = -KSKl +1 E a+1 a 

S(E\ +1 ) = -KlKl +1 E\ +l 
S{K a ) = K a , 

and extended to all of U q [gl(m\n)] by: 

S(XY) = {-) lx]lY] S(Y)S(X), 

for homogeneous X, Y G J7 g [gZ(m|n)]. Again, immediately S(K^) = K a , and 
thus 5(Id) = Id, as expected. 

S is perhaps better expressed in terms of the notation introduced for the 
coproduct A. We have, for simple generators where \a— b\ = 1: 

S(E\) = -K^ S "K^E a b . (29) 

This result is also valid for the case a = b (where — sign(a — a) = 0), and the 
formula degenerates to S(E a a ) = —K^K^E a a = —E a a , which is equivalent to 
S(K a ) = K a . Furthermore, a direct inductive prooffl shows that ( p9| ) generalises 
to the case of nonsimple generators E a t, so that is, in a sense, the most 
general expression of S. 

6 An example proof, albeit for a different definition of S, is provided in Lemma 3 of |27|. 
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3 Highest weight U q [gl(m\n)} representations 7T\ 

3.1 Introduction 

The construction of highest weight representations for U q [gl{m\n)\ involves ini- 
tially postulating a highest weight vector, which we shall call |1). The action of 
the U q \gl(m\n)\ Cartan generators on |1) is that of scalar multipliers; the details 
of these multiplications are encoded in the weight of 1 1) . Thus, if we are dealing 
with a representation labeled: 

m+n 

A = (Al, . . . , A m | A m +i, . . . , Xm+n) = Ai£i, (30) 

i=l 

where the Si are gl(m\n) fundamental weights (see § [2.3| ), then we intend |1) to 
have weight A, that is the action of K a on |1) is:[] 

/if tt .|l)=7rA(if.)-|l)^«Ml>. (31) 

Substituting q N for q, we immediately have that: 

K?-\l)=q? x -\l). (32) 

We implement the notion that |1) is a highest weight vector by declaring that 
it be annihilated by the actions of all raising generators: 

#V|l)4o, a<b. (33) 

The module V\ is then defined by the action of all possible products of the 
U q [gl(m\n)] lowering generators on |1). We may determine a basis B\ for V\, 
with elements \i) defined by: 

\i)Ap i X il X i2 ---X ip -\l), 

where the are U q [gl(m\n)] generators, p is the number of generators in the 
product, and Pi is a normalisation constant. 

We call B\ a graded weight basis, meaning that we may assign to it (i.e. to 
V\) a grading consistent with that of U q [gl(m\n)], and a system of weights (both, 
see §||). For our specific choices of A (again, see §||), V\ is finite-dimensional. 

3.2 The Kac induced module construction 

The Kac induced module construction (KIMC) is a two-stage process which 
efficiently implements the construction of B\. 

• Firstly, we construct a basis for the so-called 'even subalgebra submod- 
ule' V® CV\; this being the module of highest weight A of the U q [gl(m\n)] 
'even subalgebra' U q [gl(m) © gl(n)], viz the algebra generated by the even 
generators of U q [gl(m\n)]. That is, V£ is defined by the action of all 
possible combinations of even lowering generators on |1), where we have 
declared that |1) is annihilated by the action of all even raising generators. 

7 For consistency of the weight notation between gl(m\n) and U q [gl(m\n)], tells us 

that A a is actually the weight of |1) in terms of the U q [gl(m\n)] 'generators' E a a . 
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• Secondly, V\ is induced from V® by the repeated action of the odd lowering 
generators on V®, subject to the proviso that V® is annihilated by the 
(unique) odd raising generator, to whit: E m m+ i ■ V® = {0}. This implies 
that E a i, ■ V? = {0} for all odd raising generators E a b, a < b. Thus, we 
construct B\ from B®. 

In this process, there is a subtlety: the resultant "Kac module" (i.e. V\) 
may not be irreducible. However, we shall choose A such that n\ is a 
so-called typical representation |28), ensuring that V\ is irreducible. 



3.3 Dimension of V\ 

For arbitrary typical highest weight U q [gl(m\n)] representations V\, we have the 
following Kac-Weyl dimension formula Jl7| : 

dim(F A ) = 2 m ™-dim(V\°), where dim(V?) = TT i^tE^R (34) 

76A + 

where Aq , po and the inner product (•, •) on the gl(m\n) fundamental weights 
are presented in §2.3. 

For the specific choice A = A = (0 m \d n ), for even positive roots 7, we have: 

(A )7 ) = i £ (A, = § E (-) [a] (Aa-A 6 ) = 0, 

a<fc,[a] = [b] a<b,[a] = [b] 

as A a = A b for [a] = [b], thus dim(V^) = 1, hence dim(V A ) = 2""\ which 
simplifies things. Details of the KIMC for this case are presented in §||. 



3.4 Matrix elements 

To construct explicit matrix elements ir\(X) for a particular U q [gl(m\n)] gener- 
ator X, the action of X on each of the basis vectors of B\ must be determined. 

Whilst the action of the generators on |1) is predefined, more generally, the 
determination of the action of X on an arbitrary vector \i) requires the rendering 
of a string of generators into a normal ordering and the application of the 'KIMC 
rules' to simplify that normally ordered expression into (a multiple of) a basis 
vector. 

Thus, we must first determine an appropriate ordering (see § [4.1[ ), and then 
describe an appropriate set of generator commutations to implement that or- 



dering (viz the PBW lemma of §4.2). By the latter, we mean that we intend 



not to use all the commutators of §2.2.5 directly. Instead, we shall use expres- 



sions taken from lemmas in [£7|, |28| for the commutations between non-Cartan 
generators. 

With these tools, we proceed to build bases and explicit matrix elements for 
our particular representations in §|[ 
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4 A normal ordering and a PBW lemma 



Finding a normal ordering for a string of U q [gl(m\n)] algebra generators involves 
the recursive use of commutation relations to rewrite the string as a sum of 
strings, with respect to some chosen (hopefully natural) ordering. Both the 
initial string and the resultant may contain initial scalar multipliers, which for 
U q [gl(m\n)\ are typically algebraic expressions in q. When we speak of the 
length of a string, we shall ignore these scalars. 

A PBW lemma describes the appropriate commutations, but we must de- 
termine an ordering ourselves. Perhaps the most natural ordering is purely by 



weight (see (|24|) in §2.3), but there are reasons for choosing other orderings. 



4.1 A normal ordering for U q [gl(m\n)] 

We begin with the convention that if generators G\ and Gi are ordered, viz 
G\ ^ G2, then the string G\Gi is ordered. With this, the ordering we choose 
is based on the following principles: 

1. Our string will often be regarded as (right) acting on the highest weight 
vector |1) 6 V\, and the KIMC directs us to first build an even subalgebra 



submodule V° (see §3^2 and §5^) based on this |1), i.e. to define basis 
vectors of V® in terms of the right actions of strings of even lowering 
generators on |1). Thus, we require even generators to be greater than 
odd generators, i.e. within normally-ordered strings, even generators lie 
to the right of odd ones. 

2. Within the even generators, |1) is always annihilated by the (right) ac- 
tion of raising generators, so these must be greatest, i.e. rightmost. By 
symmetry, we then demand that the least amongst the even generators 
are the lowering generators, so the Cartan must be lie between the even 
lowering and the even raising. 



3. Within the odd generators, the (right) action of the raising generators 



always annihilates any vectors from V®, so the odd lowering must be 



lesser than the odd raising. 

4. Within the five equivalence classes created by these considerations, non- 
Cartan generators are ordered by increasing weight. Doing this ensures 
that squares of odd generators can be systematically identified and anni- 
hilated; it also facilitates a systematic way of defining basis vectors for V\ 
(see §|). 

Furthermore, (powers of) Cartan generators are ordered by index, viz 
K~f ^ Kfr if a ^ b. Doing this ensures that powers of the same generator 
may be combined. 
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We call this ordering "OL < OR < EL < C < ER" . It differs slightly from 
that (implicitly) described in |§ pl240], viz OR < OL < EL < C < ER. To 
implement it, we say that distinct weights 71 and 72 are ordered (viz 71 -< 72) 
if the first nonzero component of 71 — 72 is positive. Then, say that we are 
comparing generators G\ and G2; where 7, is the weight of Gi (see (^4|) in 



§2.3); Li is the 'lifting' of Gi, being —1, or +1 if Gi is a lowering, Cartan or 
raising generator, respectively; and, if Gi is Cartan, then let a, be its index (the 
exponent is unimportant). Then: 

A V (B A (G V (£> A {(E A F) V G)))) Gi ^ G 2 , 



where: 








A 


is 


[Gi] > [G 2 ] 


Gi is odd and G 2 even 


D 


is 


[Gi] = {G 2 } 


both odd or both even 


C 


is 


L x < L 2 


ordered liftings 


D 


is 


L\ = L2 


same liftings 


E 


is 


Li =0 


Gi is Cartan 


F 


is 


ai ^ a 2 


(implicitly) both Cartan and ordered 


G 


is 


71 ^ 72 


ordered within (implicitly non-Cartan) class 



To illustrate the ordering, for J7 g [pi(3|l)], with reference to Table]!], we have: 
i? 4 3 < E^2 < E \ < E 4 < E 2 4 < E^4 < 

V v ' " v ' 

Odd Lowering Odd Raising 

E 3 2 < E 2 1 < E 3 1 < iff 1 < K^ 2 < K^ 3 < E 1 3 < E\ < E 2 3 . 



Even Lowering Cartan Even Raising 

This ordering ensures that in the KIMC, the action of a normally-ordered 
string of generators on a highest weight vector |1) may be evaluated by the 
following ordered steps. 

1. If there are any terminal even raising generators, then the string evaluates 
to 0. 

2. If there are terminal Cartan generators, then these may be replaced by 
their known scalar actions on |1), and the string is reduced in length. 

3. Next, the action of any even lowering generators is considered. In the 
general situation, these map |1) to another basis vector of V®. For our 
particular modules V\, where A = (0 m |d n ), as V® is one dimensional, |1) 
is in fact annihilated by the even lowering generators, so if there are any 
even raising generators, then the string evaluates to 0. 

4. Next, the odd raising generators annihilate any vectors of V^, so if any 
are present, then the string evaluates to 0. 

5. Lastly, when the string is reduced to (a scalar multiple of) the action 
of some unrepeated odd lowering generators on a V® basis vector, that 
residual string may be identified as (a scalar multiple of) a particular V\ 
basis vector. 
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We mention that although our PBW lemma provides us with the means 
to normally order generator strings, the normal ordering is a computationally 
expensive process. Firstly, each exchange may generate up to two extra terms in 
a sum^j so there is a geometric increase in the number of terms with exchanges. 
Secondly, implementation of the exchanges is really a sorting procedure, but 
we have not been able to implement an efficient algorithm - we in fact use 
the dumbest possible opportunistic exchange. This failure is partly due to the 
complexities in developing a sorting algorithm in the presence of the continual 
creation of extra terms. 

Thus, the process to normally order generator strings requires time and 
storage which both of which grow at least exponentially with string length. 
Using Mathematica, we currently get into serious trouble beyond length 8. 



4.2 Commutations implementing the normal ordering 



To implement the normal ordering described in §4.1, we describe here a set 
of generator-exchanging commutations. The material originates in p8fl ; we 
have modified the results a little in light of ([l5]), rearranged many things, and 
corrected several minor mistakes. In what follows, we intend distinct abstract 
indices to represent different concrete indices. 



4.2.1 A PBW commutator lemma 

The following result contains some corrections to the original p8| . In it, we use 
the notation presented in (|Hj). 

8 Although the exchanges can in fact add up to two extra terms, in practice they add only 
one extra term, but they can be 'sum-neutral' or even subtract a term. 
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Lemma 1 

We have the following commutations. Firstly, (|7^) generalises to the case of 
nonsimple generators mw, viz: 



[E a b , E b a ] = A a {K a K b - K a K b ) all a, b. (35) 
Secondly, where there are three distinct indices, we have: 

{(a) K c K b E\ c<b<a ~) 

(b) E a b K a K c c<a<b I , , 

(c) E\K c K a b<a<c (36) 
(d) K b K c E a b a<b<c J 

[E c a , E c b ] = [E a c , E b c ] = a<c<borb<c<a (37) 

E\E\ = J W ^Z"\ 8<t<C L^ fl (38) 

\ (b) {-) lE a] c<a<b f V ' 

{[E b ] 1 

(a) (~\Ji a<b<c \q c E b c E a c . (39) 
(b) (-) lE < ] c<a<b f 

Thirdly, we describe the situation where there are no common indices, and we 
have a < b and c < d. Let S(x,y) denote the set of integers {x, x + 1, . . . ,y}. 
Then, if S(a, b) and S(c, d) are either disjoint or one is wholly contained within 

the other, viz a < c < d < b, a < b < c < d, c < a < b < d or c < d < a < b, we 
have a total of 16 cases: 

[E\, E c d ] = [E\, E d L ] = [E b a , E c d ] = [E b a , E d c ] = 0. (40) 

More interestingly, if there is some other overlap between the sets S(a, b) and 
S{c, d), viz a<c<b<dorc<a<d<b, then we have the 8 cases: 

? b Tpd 1 _ / ( a ) + A b a<c<b<d\^ d ^ 



r P a pd] _ / (a) -A b K b K c E\E d b a<c<b<d\ 

[ cJ " 1(b) +A d E d b E\K a K d c<a<d<b} 

r p0 pel _ f (a) -A c E b d E c a K c K b a<c<b<d \ , . 

L a ' dJ 1(b) +A a K d K a E\E b d c<a<d<bj- 
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Rearranging the indices in Lemma [j] gives us the following simplified results. 
• The entirety of (pq) may be summarised by: 



[E a Cl E c 



K c K b E a b c<b<a 

E a b K a ~K c c<a<b 

E a b K c K a b<a<c 

K b K c E a b a<b<c. 



• The entirety of (|37|) to (|39|) may be summarised by: 

E a c E b c = nE b c E a c and E c a E% = nE c b E c a , 



where: 



if z(a, b,c) = c 
else, 



and z(a, b, c) picks out the middle element of {a, b, c}. (The 1 factor follows 
as [£ a c ][£* c ] = for c between a and b.) 

• The entirety of (fhJ) to m4) may be summarised by: 



[E a b ,E c 



+A b E a dE c b 
—lS.dE a dE° b 
+A a E c b E a d 
—A c E^ b E a d 
-A b K b K d E a _dE c b 
+A c E%E a d K a K c 
—A c E a d E c b K c K a 
+A b K d K b E c b E a d 




a < c < b < d 
c < a < d < b 
b < d < a < c 
d < b < c < a 
a < d < b < c 
d < a < c < b 
b < c < a < d 
c < b < d < a 
a ^ b ^ c ^ d else. 



From these, we deduce the following rules for exchanges: 
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Prom (JT2J), replace Kf Kg with Kg . 



If also a = b, then replace it with K^ +N . 

If also A/ + iV = then replace it with Id. 



Prom fll5|), replace K™E b c with 9 a 



E b c Kg , and 



replace £ b c iC with qa N(K ~^ ] Kg E b c 



From fll9|), replace E\E b a with (-)^ bi E b a E\ + A a (K a K b - K a 



From fl2l|), replace E a b E a b with if [£7° 6 ] = 1. 
Replace E a c E c b with: 



(-) 
(-) 
(-) 



^E c b E a 
] E c b E a 
] E c b E a 



E\ + qc b E%E a c 

KMbE a b 
E a b K a ~K c 
E a b K c K a 
K b ~K c E a b 



Replace E c b E a c with: 



q c a {E a c E c b 



(^c^b 

(E a c E c b 



E a b )_ 

K c K b E\) 

E\K c K a ) 
K b K c E\) 



a < c <b, b < c < a 
c < b < a 
c < a < b 
b < a < c 
a < b < c. 



a < c < b, b < c < a 
c < b < a 
c < a < b 
b < a < c 
a < b < c. 



Replace E c a E c b with nE c b E c ai and 

replace E a c E b c with K,E b c E a c , where: 



1 



if z(a, b, c) 
else, 



and where z(a, b, c) picks out the middle element of {a, 6, c}. 
Replace E a b E c d with (-j^^lPd^^j _|_ j' ) where: 



T = < 



+A b E a dE c b 
—lS.dE a d.E c b 
+A a E c b E a d 
—A c E^ b E a d 
-A b K b K d E a _dE c b 
+A c E c b E a d K a K c 
-A c W d E\K c K a 
+A b K d K b E c b E a d 




a < c < b < d 
c < a < d < b 
b < d < a < c 
d < b < c < a 
a < d < b < c 
d < a < c < b 
b < c < a < d 
c < b < d < a 
a ^ b ^ d else. 
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5 The U q [gl(m\n)} representations (0 m |d n ) 

Fixing m and n, in this section we describe the use of a version of the Kac 



induced module construction (KIMC, see §3.2) in the brute- force construction 
of the U q {gl(m\n)\ representation A = (6 m |o!„). 

Alternatively, we might have implemented the results presented in j22l ^3|, 
which describe the use of a Gel'fand-Tsetlin basis to explicitly construct the 
actions for essentially typical representations^ (this class includes our represen- 
tation, which is actually typical). We avoid those fine results because we wish 
our code to be more general, but we pay a price for this in the currency of 
computational expense. 

Strictly, this material applies only to generic q, that is q not a root of 
unity (in which case the representation theory changes drastically). Also, our 
representations are unitary only under some constraints on a (viz that a is real 
and either a > n — 1 or a < 1 - to, see JTT[), and so we shall implicitly select 
these. In the application of our results to the computation of link invariants 
[^| , the representation of the braid generator based on our quantum R matrix 
R m ' n again contains the variables q and a. However, it turns out that R m ' n is 
actually a valid braid generator for any q and a. 

5.1 An orthonormal basis B for V = V\ 



Recall from § 3.2 and that V = V\ is of dimension 2 mn , and may be equipped 
with a grading compatible with that of U q [gl(m\n)]. It is known that V contains 
no weight multiplicities (that is, V contains no constant weight subspaces of 
dimension greater than 1), so that a weight basis for V will contain no distinct 
vectors of the same weight, and this makes our task a little simpler. 

Here, we use a version of the KIMC to construct a weight basis B — {|«}}f =1 
for V. That is, the 2 mn basis vectors \i) are defined in terms of the actions of 
all 2 mn possible nonrepeated, ordered combinations of the U q [gl(m\n)] simple 
lowering generators (there are mn of them) on a postulated highest weight 
vector |1). This |1) is further defined to be of unit length and annihilated by 
all U q [gl(m\n)] raising generators. Each of the vectors defined in this manner 
will be orthogonal to all other such vectors, and it is a straightforward matter 
to select constants to orthonormalise them. The resulting B is thus a graded 



orthonormal weight basis for V . Using it, in §5.2, we construct matrix elements 
for the Uq[gl(m\n)] generators. In each subsection, we shall illustrate our results 
using U q [gl(3\l)}. 

9 It turns out that these results sometimes hold for other representations, when various 
limits are evaluated using L'Hopital formulae. 
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5.1.1 Details of the KIMC 



The KIMC firstly instructs us to construct a basis B° for V°, the submodule of 

V determined by the action of the U q [gl(m\n)] even subalgebra on |1). For our 
choice of A, in fact V° has dimension 1 (see § |3.3| ), hence B° = {|1}}. 

This means that |1) is annihilated not only by all raising generators, by 
also by all even lowering generators. More generally, for other representations, 
we have to work harder to construct a basis for V°; that process has a similar 
appearance to the following. 

Secondly, B is induced from B° by the actions of all possible products of 
odd lowering generators on B . Thus, we must consider the set of all possible 
products of odd lowering generators. The PBW lemma allows us to reduce 
this set to that of all possible ordered products, and the knowledge that the 
square of odd generators is zero allows us to reduce it to the set of all possible 
nonrepeated ordered products, a finite set, ensuring that B is finite. Recalling 
that U q [gl{m\n)\ has mn (simple and nonsimple) odd lowering generators: 

OLGS = {E a b :a = m+l,...,m + n,b=l,...,m}, (45) 

thus V is spanned by a set of vectors obtained by the actions of all possible 
ordered products (that is, strings of length to mn) of OLGS generators on 
|1). Indeed, this is the source of the factor 2 mn in the dimension formula ([34]). As 

V is known to have no weight multiplicities, this spanning set is itself the desired 
weight basis B. The set of ordered products of generators may be obtained from 
the power set V(OLGS), by replacing its elements with respective products - 
Mathematica is well-suited to this. 

To illustrate, for the U q [gl(3\l)] case, we have OLGS = {E 4 3 , E 4 2l E 4 ^, 
and V has the following basis B: 



/ 


1) 


\ 


/ 




\ 




2) 












3) 






p 3 ES 






4) 






PiE\ 






5) 






P$E A zE A 2 






6) 






P 6 E 4 3 E 4 i 






7) 






(3rE 4 2E 4 i 




\ 


8) 


/ 


V 


PsE 4 3E 4 2 E 4 i 


) 



where the /Jj are scaling factors, which we shall select to normalise the the 
redundant (3\ is implicitly 1. 
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5.1.2 Weights, gradings and an ordering for B 



Our B is a weight basis, in that each of the \i) is of a definite weight. If we have, 
for some i: 



(47) 



\i)=fcE a i bl E a \ 2 ---E a * bp -\l), 

for some p ^ mn, then we may define: 
p 

wt(|i))=A + ^wt(^ 6j .), 
i=i 

where wt (£?"(,) is the weight of the generator E a \> (see ( p4]) in § [2.3| ). As we 
defined |1) to be a highest weight vector, clearly we intend wt(|l)) = A. 

Further, our B is Z2 graded, in that its elements are formed from the actions 
of products of graded U Q \gl(m\n)] generators on the (zero) graded |1). The Z2 
grading of \i) (as defined in (E7n) is defined by: 



[[1>] = 0, [|i>] 



p (mod 2), 



where the latter result holds as the E aj b j are all odd. This Z2 grading on 
V is compatible with a notion of Z graded level, this being p, the number of 
factors in the product forming \i). This notion is relevant in the calculation - 
we recursively form the basis vectors in level I by the action of the OLGS on 
the basis vectors of level I — 1, for 1 = 1,..., mn. We number our vectors \i) by 
decreasing weight within increasing Z graded levels. This ordering is important 
in that it simplifies t he p rocess of identifying an arbitrary string acting on |1), 
which is required in |5.2| . 

The weights and gradings of the basis vectors for our C/ g [gi(3|l)] example 
are supplied in Table | (cf. Table @). 



i 




wt(|i)) 






[H>] 




1 


( 0, 


0, | 





) 








2 


( 0, 


0,-1 | 


a 


+ 1) 


1 




3 


( 0, 


-1, 1 


a 


+ 1) 


1 


1 


4 


(-1, 


0, 1 


a 


+ 1) 


1 




5 


( 0, 


-1,-1 1 


a 


+ 2) 







6 


(-1, 


0,-1 1 


a 


+ 2) 





2 


7 


(-1, 


-1, 1 


a 


+ 2) 







8 


(-1, 


-1,-1 1 


a 


+ 3) 


1 


3 



Table 2: Weights and gradings of the basis vectors \i) of B, for {7 g [flrZ(3|l)], 
ordered by decreasing weight within increasing Z graded levels [|«)] z - 
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5.1.3 Normalisation of B 



To investigate questions of orthogonality, we require an inner product on V. To 
whit, we introduce a basis B* = {{i\}f2\ 01 ^* (the dual of V), by: 

<i|4/^(l|.E fe % p --.E 6 %E bl 0l) (48) 

where ft is the complex conjugate of ft, and \i) is as supplied in ([f7|). Note 
that here, we explicitly intend the Hermitian conjugates E b a (see (|J)), and not 
E b a . This ensures that B* is conjugate to B, viz we have \if = and that 
conjugate generators and conjugate basis vectors remain conjugates in their 
matrix representations. We assign weights and Z2 (and also Z) gradings to the 
(i\ such that wt((t|) = wt(|i)) and [(t|] = 

Using this conjugate basis, we define an inner product on V: 

(\i),\3)) = (i\-\3) = m, (49) 

where we implicitly have (1|1) = 1. Next, where E b a is any raising generator 
(viz b < a), taking the conjugate of E h a • |1) = yields (1| • E a b = 0. Expanding 
the E a b into simple generators using (|b) shows that (1| • E C+1 C = 0, for all 
simple lowering generators E c+l c , which in turn, recursively, yields (1| • E c d = 
for all lowering generators E c d with c > d. In sum, the equivalent of (|3^) is: 

(1\-E\ = 0, a>b. (50) 

For completeness, we mention the (left) action of Cartan generators on (1|, 
obtained by conjugating (|3^): 

(l|.Jf^ = g^<l|. (51) 

More generally, the value of an inner product may be calculated by the 
following procedure. 

1. We substitute the definitions of (i\ and \j) (viz the appropriate versions 
of © and @) into ©, yielding a form ft fa (1| • Z ■ |1), where Z is a 



string of generators. For our representations tta, as V is one-dimensional, 
Z will contain no even generators. 

2. We use the PBW lemma to normally order Z . The resulting strings have 
their raising generators annihilating |1) pushed to their right hand ends 
and lowering generators annihilating (1| pushed to their left hand ends. 

3. Implementing those annihilations, and evaluating the residual (scalar) 
Cartan generator actions on |1) and (1| (viz ( |3^ ) and (|5l])), we convert 
the resulting expression to a scalar. 
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It turns out that vectors \i) and |j) with distinct weights satisfy = 0, 
viz — S'j Recall that our (0 m |d„) representations have no weight 

multiplicities, thus distinct basis vectors have distinct weights, hence our basis 
B is orthogonal. To make it orthonorma/, we must select the /3j appropriately. 
This means that for each i, we must ensure that: 



(i|i)=/3*A(l|-E 6 



■ ■ ■ E 



,E 6 



E ai bl E 



Thus, for each i, we must use the commutations of the PBW lemma to normal 
order the following string Z{. 



Zi = E° 



• • • E 



,E 



In 



E" 



■ ■ ■ E ap 



into NO(Zj), and then apply the algebra- module actions fl32|), (|33|), (|50|), ( pl| ) 
and (1 1 1) = 1 to (1 1 • NO(Zj) • |1) to yield, up to an arbitrary complex constant 

— 1 /2 

(phase factor), = ((1| • NO(Z,) • |1)) . The phase factor is unimportant; 
different choices simply lead to bases related by orthogonal transformations, and 
this will not affect our R matrices. In practice, we let the internal machinery of 
Mathematica decide on phase factors for us - a human calculator might make 
more elegant choices. 

So, at this stage, we have determined the constants /Jj such that we have 
an orthonormal basis B for V . In general, for arbitrary representations of 
Uq\gl(m\n)], these constants will be algebraic functions of the complex variable 
q. For our particular A, these functions will also contain the variable a. 

To illustrate, for the t/ g [gZ (3 1 1)] case, we have: 



ft 




[a 


" X 




ft 


= q 


[a 


»: 




ft 


= t 


[a 






ft 


= q 


[a 




[a 


ft 


= t 


[a 




[a 


ft 


= t 


[a 


q i 

_ 2 


[a 


ft 


= t 


[a 


q i 
7 


a 



} (52) 



+ 2]" 5 



where we have use the g-bracket (see (|1 
for example, 



to simplify the expressions. Thus, 



ft 



q 6 (q a -q a )(q aA 



{q-qf 



Our use of the g-bracket notation is more than cosmetic; the U q [gl{m\n)\ 
symmetries manifest themselves naturally in n\ in these patterns, and if we do 
not recognise and incorporate them into our notation, then expressions rapidly 
become unreadable, and then intractable. Below, g-brackets will appear at every 
point, and even in our R matrices. 
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As mentioned in §4.1, the normal ordering of generator strings is a computa- 
tionally expensive task. Here, the normal ordering of Zymn, typically dominates 
the computations as it demands that we process a seriously disordered string 
of length 2mn. A theoretical insight would be valuable here - for example an 
explicit formula for the normal ordering of arbitrary E b a E a i, would help speed 
the evaluation of the $ . The regularities apparent in the above example suggest 
such the existence of such a result, and (Q) may also be of use. Alternatively, a 
more efficient computation of the should be possible by the efficient reuse of 
previous calculations. 
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5.2 Matrix elements for tx = 7Ta 

Having established an orthonormal basis B for the module V = Va correspond- 
ing to the representation 7r = 7Ta, we now use it to construct matrix elements 
7r(X) for U q [gl(m\n)] generators X. 

Where /2 m ™ is the identity transformation on V, in the basis B we have the 
identity: X)»=i = ^2™™, thus: 

2 m " 2 m " 2 m " 

tt(X) = tt(X) • ^ |i)<i| = 5>(X) • • (i| = • |i)) • (i| , (53) 

i=l i—1 i—1 

and the action X ■ \i) may be computed knowing the expansion of \i) in terms 
of generator products and |1). Thus, if, as in (|47|), we have: 

\i)=0 i E a \ 1 E a * b2 --.E a » bp -\l), 

then we may compute X ■ \i) by the following process: 

1. We again use the PBW basis commutators to normally order the string 
Y = XE ai bl E a2 b2 ■■■E a " bp , denoting the result by NO(Y). 

2. We use the known actions of the raising and Cartan generators on |1) to 
reduce NO(Y") ■ |1) to an expression which is generally a sum of scalar- 
multiplied, normally ordered products of (odd) lowering generators acting 
on 

3. Identifying the terms in the resulting products as scalar multiples of var- 
ious \j), we obtain the result. To whit, if NO(Y") contains a term of the 
form: 

T ■ \1) = 8E Cl dl E c ' 2 d2 ■ ■ ■ E Cr dr ■ |1>, 
for some scalar 9, and r odd lowering generators E°i , and we know that: 

\k)=l3 k E c ' dl E C2 d2 ---E^ dr .|1>, 
then we may replace T ■ |1) with Q(3 k \ k). 

Repeating this procedure for all 2 mn basis vectors and substituting the 
results into ([53]) yields the required matrix element 7r(X). Again, as the ordering 
chosen for the basis vectors |z) is compatible with the ordering used in the PBW 
lemma, this process is robust. 

We now divide the construction of matrix elements into two phases. Firstly, 
in a direct implementation of the above, we build matrix elements for the Cartan 
and simple raising generators. We illustrate this for the t7 9 [(7/(3|l)] case, for the 
generators K4 (in |5.2.1 ) and E 3 ^ (in §5.2.2). 

Having done that, in § [5.2.5 , we describe the construction of the remaining 
matrix elements, as they may be efficiently computed in terms of those for the 
simple lowering generators. 

Further illustrations are provided in my PhD thesis B, using the C/ g [<?/(2|l)] 
case, although those results are somewhat less formally explained. 
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5.2.1 Matrix elements ir\(K4,) for the C/ g [g^(3|l)] case 

Firstly, we must normal order a list of 8 generator strings, cf. (fhf). We obtain: 



/ Id 



\ 



E\ 

T7l4 jp4 

Ei iHi 2 

TP 4 TP 4 
Hi 'i£j l 

TP 4 TP 4 
Hi 2£j 1 

r,4 _ t-.4_ 



/ 



(15) 



\ iS 3_B 2-E i / 



Id 

qE 4 3 

4 

2 

1 

1 



qE 4 - 



qE 4 - 
q 2 E 4 3 E 4 2 



q 2 E 4 3 E 4 , 
q 2 E 4 2 E 4 1 
\ q 3 E 4 3 E 4 2 E 4 1 J 



(54) 



The action K±- 11) is known explicitly from 
Thus, we have: 



2h, that is, we have J^" 4 -|l> = <f*|l). 



/ 



K 4 



\ 



!) \ 

2> 
3) 
4) 
5) 
6) 
7} 

8> ) 



(P) _Q 

u 3 



1 

WiE 4 i 
i 

2 



qfaE 4 ! 
q 2 feE 4 3 E 4 2 
q 2 06E 4 3 E 4 ! 
q 2 PlE 4 2 E 4 1 
\ g 3 /3g_E 4 3_E 4 2-E 4 i 



1) S g° 



/ 


11} \ 




|2> 




|3> 




|4> 


<f 


|5> 


-2 

<7 


|6) 




|7> 


U 3 


|8> / 



Installing this information into (|53T), we discover 

7T(# 4 ) 



9 



|l)(l|+g(|2)(2| + |3)(3| + |4)(4|) 

+g 2 (|5)(5| + |6}(6| + |7)(7|) + g 3 |8)<8| 



Note that this process doesn't actually require explicit knowledge of the fli 
presented in (52). This is a feature of the evaluation of matrix elements of 
Cartan generators; the more usual situation appears in the next example (i.e. 
§ 5.2.2p . Replacing with the elementary matrix e*- <G we have: 



that is: 



n(K 4 



(el- 




(el 
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4 + 


et)- 




(4 
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7/ 
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5.2.2 Matrix elements 7r A (-B 3 4) for the U q [gl(3\l)] case 



This is a more interesting case than that of K4 . 
of the list of 8 generator strings yields: 



/Id 



E 4 3 



E\ 



\ 



TP 4 TP 4 
Hi 3 Hi 2 

TP 4 TP 4 
Hj 3 H, l 



? 4 2 

P 4 . cvi 



J? 4 IP 4 
Hj 2H1 1 



\E\E\E\J 



(E Z A 

-E\E 3 4 + A(K 4 Ks - K 4 K 3 ) 
— _B 4 2-E 3 4 + E" '2K4K3 
-_E\E 3 j_+ E z 1 K i K 1 _ 
AE 4 2 (qK 4 K 3 - qK 4 K 3 ) + E 4 3 (E 



This time, the normal ordering 
\ 



4 ( TP 4 jp3 
2Fj 4 

E 4 3 {^E 4 \E 3 4 



AE\(qK 4 K 3 -qK 4 K 3 _ 
_E 4 2-E 4 i-E 3 4 + qE 4 1E 3 2K 4 K 3 
-E 4 3 (E 4 2 E 4 1 E 3 4 + qE 4 1 E 3 2K 4 K 3 



E 3 2 K 4 K 3 ) 
E\KiK 3 ] 



E 2E \K 4 K 3 



E 4 2 E 3 1 K i K 3 ) 



\ 



+AE 4 2 E 4 1 (q 2 K 4 K 3 -q 2 K 4 K 3 ) ) 



where we have again written A = q — q and A = A -1 (cf. (|Io|)), and done a 
little judicious factoring to improve readability. Next, we know that any terms 
which end with raising generators (e.g. E 3 4 ) will annihilate |1). Further, any 
term which contains an even lowering generator immediately to the left of a 
terminal string of Cartan generators^ will annihilate |1), as the action of the 
Cartan generators is purely scalar. Omitting all such terms, we may thus write: 



E\- 



Id \ 

eS 

E 4 2 

E 4 1 

TP 4 TP 4 
Hi 3 Hi 2 

E 4 3 E 4 , 
E 4 2 E 4 1 
\ E 4 3 E 4 2 E 4 1 ) 



( 



|1> = 



A(K 4 K 3 - K 4 K 3 



\ 

3) 





AE 4 2 (qK 4 K 3 -qK 4 K 3 
AE 4 1 (qK 4 K 3 ~qK 4 K 3 







11) 



\ AE 4 2 E 4 1 (q 2 K 4 K 3 -q 2 K 4 K 3 ) J 



Again, the action ■ |1) is known explicitly from (p2|), that is, we have 
Kf ■ |1) = |1) and, as above, that K± ■ |1) = q Ta Thus, we have: 



E\- 



Id \ 



E 4 3 



TP 4 TP 4 

Hi 3 Hj 2 

E 4 3 E 4 4 
E 4 2 E 4 4 

4. c4. 



\ E\E\E\ ) 



I 



|1> = 



\ 






[a + l] q E 4 2 
[a + 1]E 4 1 



V [a + 2] q E 4 2 E\ J 



|1>- 



(55) 



10 For example, the last product contains E 4 3 E 4 2E 3 4 K 4 K 3 , which contains the even low- 
ering generator E 3 \ immediately to the left of a terminal string of Cartan generators. 
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This time explicitly installing information of the Pi from (bin, we have: 



E\ 



|2> 
13} 
14} 
1 5) 
1 6) 
1 7) 
VIS)/ 



(40) 
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faE 4 - 



f3^E 4 3E 4 2 
^ 4 3 £ 4 i 
l3rE 4 2E 4 i 
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ID© 



(0) 



foE 

\f3%E 3E 2E 1/ 

#»[«], 
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^ a [a + l] q \3) 
/3 6 /3 4 [a + l]j4> 


\/3 8 /3 7 [a + 2]j7>/ 



/? 6 [a + l] ? £ 4 a 

V/3 8 [a + 2] ij .E 4 2 .E 4 1 

0\ 

1) 





°\ 

IX) 



ll> 
II) 


IX)/ 



|3> 

|4) 


V[« + 2]| \7)J 



[a+l]| 
[a + l]| 



(56) 



Note that there is a subtle point in the second application of (|46|), in that 
we are using it implicitly: this point wasn't so clear when we computed tt(K4). 
Whilst this is simple enough for a human to perform, computer programs require 
explicit algorithms. In practice, what this means is that we must invert (ff6|) to 
provide a list of transformation rules for strings of odd lowering generators acting 
on |1) that allow us to recover the \i). These rules must be carefully coded, and 
applied in reverse order to (|f6|), to ensure that we collect the longest strings 
first. Again, the ordering of the \i) facilitates this. To illustrate, we want to 
apply, in order, the following rules: 



( E 14 B 14 P 4 


|1> ► 


- 08 


18} ) 




|1> - 


+ 07 


17} 


» 3 fi 1 ' 


|1> - 


- 06 


|6) 


S 14 F 4 


|1> - 


- 05 


15} 


E\- 


|1> - 


+ 04 


14} 




|1> - 


- 03 


13} 


{ E\- 


|1> - 


02 


12} J 



Installing the information in (56) into (p3l), we discover, again substituting 



e) for \i)(j\ 
*{E\) 



|1)(2| 



l]J|3)(5| + [a + l]||4)(6| + [a + 2 ]J|7}(8| 



]J4 
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and interpreting e* as an elementary matrix: 



/ 
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[a + 2)1 
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5.2.3 Matrix elements for the remaining generators 

Having found the matrix elements for the simple raising generators 7r(_E a b), we 
may immediately write down the matrix elements for the corresponding simple 
lowering generators ir(E b a ), as these are necessarily the transposes of ir(E a b). 
To illustrate, for the [/ q [<7/(3|l)] case, we have: 



2]|e?. 



Beyond this, we may construct matrix elements for the nonsimple generators 
from those for the simple ones by (recursively) applying tt to (|^a), viz: 

ir(E\) = n(E a c )7T(E c b ) - q?n{E c b )n{E a c ), 



although this may not yield particularly useful results. Note that we can in 
fact write down the matrix elements for any (simple and nonsimple) lowering 
generators directly from their corresponding raising generators, as we have, for 
a > b: ir(E a i>) = Tr(E b a ) Tq , where the 'q transpose' T q indicates the combination 
of the transpose and the mapping q i— > q. This follows from inspection of the way 
that (||a) depends on q; indeed it is trivial for simple generators. Illustrations 
are visible below (and in Appendix B.l); the reader should note that the q 
bracket is invariant under q i— > g. 

For completeness, we list the matrix elements of the C/,j[(^(3|l)] simple gen- 
erators (the matrix elements of all the generators are listed in Appendix B.l): 



(a) 


7r(Xl) 


= e\ + e\ + 


(b) 


tt(K 2 ) 


= e\ + el + 


(c) 




= e\ + qel 4 


(d) 


HK A ) 


= q a e\ +q a 


(e) 


%{E X 2 ) 

7T(E\) 

n(E 2 3 ) 


= -4 - 4 


(f) 


- -p 4 - p 6 

— e 3 e 5 


(g) 


— e 3 e 7 


00 


tt(E\) 


= -4 - 4 


(i) 


ir(E\) 


= [°\\4 + [ 


(j) 


<E\) 


= N|ef + [ 



,3 _: 77^4 



1^3 
3 
3 



qe\ + e$ - 
-ei + qel 

4 
4 



e^ + qel + qel 



■ q4 + 9 e ? + 14 

" ' qe 7 7 + qel 



4 + 1 e 8 



\+q a+1 {4 + 4 + ei) + q a+2 {e 



^4 



l]^4 + [a- 



e 7 7 ) + q a+3 e 8 8 



•2]|ef 



} (57) 
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6 The submodules V k C V ® V 



We now turn our attention to the tensor product module V ®V . Where V 
has a basis B = {|i)}| =1 , the 2 2mn dimensional V <E> V has a natural basis 
{\i) ® \j)}ij = i, which inherits a weight system and a grading from V: 

wt(|i)®|i» = wt(|i))+wt(|i)) (58) 

[|i>®|j)] = [|i>] + [|i>] (mod 2). 

To build an R matrix acting on V®V , we will use an alternative, orthonormal 
weight basis 58 for ® V, which corresponds to the (known) decomposition of 
V ®V into irreducible U q [gl(m\n)] submodules. The basis vectors of 23 are 
expressed as linear combinations <8> \j)), where the coefficients 8ij are in 

general algebraic expressions in q and a. Before proceeding, we introduce some 
machinery for dealing with the tensor products. 

6.1 Tensor product representation tools 

We first introduce the action of U q [gl(m\n)) <E> U g [gl(m\n)] on V ® V. Where Y 
is an homogeneous U q [gl(m\n)] element, we define: 

(X ® Y) ■ ® |j» 4 (-)MM (x • |») ® y ■ (59) 

and extend by linearity to all of U q [gl(m\n)] ® U q [gl(m\n)]; note that we have 
written [i] = for readability. This action is compatible with the grading. 

Next, we define a basis {{i\ ® (il}f,7=i for ® V)* (the dual to V ®V) by 
dualising each of the elements of the basis {\i) ® ™=i : 

<*l ® 0*1 = C-) Wb1 (|i> ® U» t » (60) 
where we intend ® Ij'}) 1 = |i) t ®|j) t . The conjugate is extended by linearity: 
(A \i) ® |j) + B ® ||»t 4 A*(|«) ® Ij}) 1 + B*(\k) ® ||))t, 

for scalars A and £?, where A* is the complex conjugate of A. 

The multiplication operations between the dual bases are given by: 

(ii)®ii»-«*i®<ii) = n m] i\m®\3)(i\) (6i) 

«i|®0'IH|fc>®l*» = (-) im SlS{, (62) 
analagously to (|59|). 

Thus, we have the natural inner product on the basis {\i) ® b')}f ,-=i^ 

® li) , |fc> ® 10) - (10 ® li)) f • (1*0 ® 10), (63) 
which behaves as expected: 

(10 ® li) , I*) ® 10) ® (-) [il[j1 (W®01)-(|fc)®IO) 

= (-) m {-) im ({i\k}®{j\l)) 
Lastly, we will often use the shorthand |i® j) = \i) ® |j) and (i®j'| = (i| ® 



35 



6.2 The orthogonal decomposition of V <S> V 

For our modules V = V\, the orthogonal decomposition of V <8> V is known, and 
contains no multiplicities JO], |l5). To describe it, we introduce a little notation. 
What follows is strictly true only for m ^ n, but the natural isomorphism 
between U q [gl{m\n)\ and U q [gl(n\m)] shows that we need not consider the other 
case n m. 

For any ^ N ^ mn, let 7 be a nonincreasing sequence of nonnegative 
integers^] 71, 72, . . . , j r , where the 7^ satisfy J2k=i 7* = N. We then define a 
Young diagram [7] = [71, 72, ... , j r ], to be allowable if it has at most n columns 
and m rows, viz r ^ to and, for each k, 7& $J n. To each such allowable diagram, 
we associate a gl(m\n) weight A 7 : 

A 7 = (6 m _ r ,-7 r , ...,-72, -7l I »V,(r- l) 7ri _ 7r , . . . , i 7l - 72 , 6„_ 7l ). (64) 

Then, for our specific representations with A = (0 m \o-n), modulo the comments 
on limiting a in §|^ (viz a must be real and either a>l-nora<l - to), we 
have the following irreducible decomposition of Va: 

[7] 

where the sum is over all possible allowable Young diagrams, and V® +x is a 
U q [gl(m) © gl{n)\ module of highest weight A + A 7 and Z graded level N. More 
interestingly, we also have the following decomposition: 

V A ®V A =0F 2 a+a t , (65) 

[7] 

where again the sum is over all possible allowable Young diagrams, and V2A+.V, 
is a U q [gl(m\nj\ module of highest weight 2A + A 7 and Z graded level N. 

And now, we introduce an abuse of notation. Instead of the explicit (|65|), 
we shall often write: 

V®V = @V k , (66) 

fe 

where the submodule V k has highest weight A^, and the generic index k runs 
over some appropriate index set. In the special case n = 1, the decomposition 
of ( |65| ) becomes: 

m 

V®V = @V k , (67) 

k=0 

where V k has highest weight X k = {0 m -k, (— 1)*. | 2a + fc), thus the submodules 
Vo, Vi, . . . , V m are ordered by increasing Z graded level (which, in this case, is 
in fact, k). We shall be using (|67|) to illustrate specific examples later. 



L We apologise for overloading 7 - this one is a sequence, not a root. 
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In the following subsections we describe the construction of orthonormal 
weight bases 03 fc = )}f^ Vk) for each V k , where span(U fc 03 fc ) = V ®V. 
Fixing k, each j^f) is a linear combination of terms of the form 9ij As 
we may safely define the duals (^>f | = l^f) (as expansion of the conjugate is 
possible using (BQ)), it is thus meaningful to orthonormalise the j^f). 

Firstly, in § |6.3| , we determine a highest weight vecto r K ^i) for which 



we use as the starting point for the KIMC. Next, in £6.4 we build a basis 
03°, for the even subalgebra submodule V fc °. In this distinct from that 

of V° in Vj^' may not be one-dimensional (although it contains no weight 
multiplicities), so this construction is nontrivial, although it turns out to be 



quite straightforward. Lastly, in §^5| and §6.q , we construct 03fc by the actions 
of the 2 rnn possible combinations of ordered nonrepeated lowering generators on 
03°,. A subtlety in this case is that Vk in general contains weight multiplicities, 
so that we must employ a Gram-Schmidt process to orthonormalise it. 

6.3 A highest weight vector hj/*) for Vk 

We begin the construction of 03 & with the deduction of a highest weight vector 
of weight Afe. As we know the weights of the \i), using (|58|), we may 
immediately write down: 

|*i> = E i A'N®i>, (68) 

where the sum is only over i, j such that wt(h&f)) = = wt(\i)) + wt(\j)) is 
satisfied (i.e. we don't know in advance how many terms there are in the sum). 
The coefficients 9ij are scalar expressions in q and a, which we shall determine 
by the following. 

1. We demand that be annihilated by the actions of (the coproducts 
of) all raising generators. (Actually, it is sufficient that it be annihilated 
by all simple raising generators.) As dim(VJ?) 1 in general, j^i) may 
not be annihilated by the actions of even lowering generators. 

2. We demand that j^i) be normalised. Examination of ( |68| ) shows that, 

using (|6|), the magnitude of is (X)ij(-) [<1[jl #*j%) 2 > where the sum 
is again only over appropriate As the #y are only determined up to 
arbitrary phase factors, we may replace OhOij with 0?-, hence we demand: 



£ y (-r l ^, = i- (69) 

It turns out that these considerations always yield exactly enough constraints 
to uniquely determine j^x) (well, up to an unimportant phase factor). 

Before proceeding, we observe that j^i), the highest weight vector of the 
first module Vo, is necessarily |l<g) 1). This provides a check on the methods used 
to determine the other More generally, we must set up a linear system to 

determine the coefficients 9ij, and again, Mathematica is well-suited to this. 
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6.3.1 Illustration: |#?) for U q [gl(3\l)} 

We now use Hp?) for C/ g [gZ(3|l)] to illustrate the entire process. Using the 
notation of (|37]), we have 4 submodules their highest weights, dimensions 
(obtained from (|34|)) and suitable l^f) are presented in Table |[ Note that for 

n = 1, for A = (6 m |a), in fact ( p4| ) degenerates to dim(V°) = ( ™J . 



dim 



k 








A fe 






v? 


v k 





( 


o, 


o, 





2a 


) 


1 


8 


1 


( 


0, 


o, 


-1 


2a- 


-1) 


3 


24 


2 


( 


o, 


-1, 


-1 


2a- 


-2) 


3 


24 


3 


(- 


-1, 


-1, 


-1 


2a- 


-3) 


1 


8 



l*f) 



1711 

9 n |1®2)4 

0i 5 |1®5) 
-0 23 12(8)3) 



|1®1) 



021 |2 

051 

032 



6»is 1 1(8)8) + 6»si 

+0 27 1 2 <g> 7) + 6*72 
+036 1 3(g) 6) + 6*63 
+6> 45 |4®5) + 6*54 



>1> 

5®1) 

3(8)2) 



88)1) 
78)2) 
68)3) 
58)4) 



Table 3: The tensor product submodules for the J7 g [gi(3|l)] case. 
From Table |, for |*f), we must determine 4 coefficients: 

I*?) = 015 |1®5> + 061 15®1> + 6»23 |28>3) + 6*32 13®2> . (70) 

The U q \gl (3 1 1)] simple raising generator set is SRGS = {E 1 2 , E 2 3 , E 3 ^}, and 
from (j2qb), the coproducts are: 



A^a) = {E\ ® K\ Ki ) + {KIK\ 



A(E\) 



(E 2 



E\) 
E\), 



(71) 



<KjKi) 

A(E 3 4 ) = {E 3 A ®KlKj) 

so the actions that we want are: 

A(E\) \ ( 

A{E\) 
A(E 3 4 ) 

To evaluate these products, we take the known matrix elements of the un- 
derlying representation J5^), and substitute these into the evaluations of the 
coproducts (|7~i|). Thus, for the example A(E 3 4) ■ |\l/?), we find: 




A(E 3 4 )-\^) 



E 3 



(E 3 4 (8 Kl K'l ) + (JT| K\ w n a) 
(015 |1> ® |5) + 051 |5) ® |1) + 23 |2) ® |3) + 3 2 |3) ® |2)) . 



38 



To illustrate the multiplication: 



E 4 <g> iff Kl 



ID)© 



(|5)®|1})^^4-|5) 



|3)®|1> 



At this point, to save space, we introduce a little more notation, which 
eliminates the q brackets altogether: 



1 

Ci. 



= ([ia+j] q ) z , where z G {±, 1} 



Ja+j 



-ia+j 



(72) 
(73) 



where i,j ^ 0. In these expressions, if i = 1 we shall simply omit it, that is, we 
intend: Cj = Cij = q a+ i + q a+ i and A' = Af^ = ([a + j] q ) z ■ Occasionally, we 

will write A*j = A~ 7 - and C ltj = (Cjj) -1 . 
Using this notation, we have: 



A(£ 3 4 ) 



tf ?) = (j£ 3 4 ® ^1^4 ) + {Kl Kl ® E\) 

(015 |1) ® |5) + 051 |5> ® |1) + 23 |2> ® |3) + 32 |3> ® |2)) 



(05i g f A? 



*32<? ; 



(|1>®|3» 
-*4)(|3>®|1» 



and, altogether, for the three SRGS generators, we have: 

/ o 

3 4) 



,A(E 3 - 



(-0239^- 032^) (|2)®|2)) 
15^4 +0239 f + ^|) |3)) 

y \+ [e^Ai - 032^+^4) (|3) ® y 



As each component of the RHS must be zero, and the \i) (8 |j) are linearly 
independent, we thus obtain a net 3 linear constraints on the 0^ from this set: 



-q 2 t>23 







M o 2 023 = 



015 + Q" 
5 f 4^51 

better written in matrix form as: 



q* + 5^032 



V 





gf^J 




(7 2 + 2 





-q 



-q 2 








/ 015 \ 

051 
023 
V 032 J 



(74) 



Note that we have left the system exactly as supplied by the raw action of the 
raising generators. A human calculator might delete superfluous signs, and per- 
haps do some factorisation, but our M ATHEMATICA code would require explicit 
instructions for this finicky and unnecessary work, so we omit it. 
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Thus, application of the requirement that be annihilated by the simple 
raising generators yields a linear system of 3 equations in 4 variables. A final 
constraint to completely determine the variables is now obtained by requiring 
that be normalised. 

In our example, using ( |69[ ) and (|7C|), and recalling the gradings of the basis 
vectors (in Table |2|), we thus have the nonlinear constraint: 

0? B + 06i-*a3-0aa = l. (75) 

Combining the information in (frij) and ((75|), we are able to solve for the un- 
knowns Oij, uniquely up to the usual overall phase factor. In practice, we may 
actually avoid the use of (fT5|), by first feeding (Q) to the Mathematica equa- 
tion solver, which returns us an answer with a free parameter (the first unknown 
#y)-P3 Setting that free parameter to unity, we obtain a suitable unnormalised 
h&i), which we may immediately normalise. To illustrate, we find: 

i i , i ,i 

TjQ+l A 2 n a+1 A 2 773 4 2 n2 A 2 

I*?) = |1®5) + |5®1> - |2®3> + |3®2> . 

r* 2 4 2 c 2 A 2 c 2 A 2 c 2 A 2 

°1 A 2,l °1 /1 2,1 °1 /1 2,1 °1 A 2,l 

To complete the results, we have obtained the highest weight vectors for each 
of the 4 submodules for the t/q[gZ(3|l)] case. These are, after a little factoring: 

= |1®1) 

=c\ [ft |2®1) -gf [1(82)) 

|*?) = CfA~li (q a+1 Al |1®5) +g Q+ M| |5®1> -gi^? |2®3) + g*Aj |3®2> 

|*?) = Cfc|3| 3 x 

/- g^ +3 A| jl®8> + g^f +3 a| |8®1) + gf +2 A| |2®7) - gT+ 2 AJ |7®2)\ 
\ - gt +1 Af |3®6) + gf + MJ |6®3) + gt^J |4®5) -g*A| |5®4) /' 

Observe the presence of q graded symmetric combinations of |i<2)j) and |j(8>i) 
in these expressions, viz patterns of the form q x \i®j) ±q x \ j®i). This feature 
appears repeatedly throughout the bases for the Vk- 

So, at this stage, we have described how to construct normalised highest 
weight vectors h&i) for each Vk- An interesting outstanding point about our 
process is that the two demands that be annihilated by the m + n — 1 simple 
raising generators and that it be normalised, yield exactly enough constraints 
to determine it uniquely (up to a phase factor) . The reason for this balance lies 
buried in the combinatorics of the ways the weights of the underlying module 
V can be added to yield the weight Afe. 

12 In fact, it returns us two answers, differing by a (phase) factor of —1. We judiciously 
choose to ignore the second one. 



40 



6.4 A basis 03° for F fc ° 

Having determined we now apply the first stage of the KIMC to construct 

the basis 03^. That is, we construct basis vectors \ by the repeated action of 
the m + n — 2 even simple lowering generators (a set which we call ESLGSf^ 
on the known h^i). In our case, as V® contains no weight multiplicities, vectors 
of distinct weights created in this way will naturally be orthogonal. At each 
stage, we must also check to see if newly minted vectors are scalar multiples of 
previously found ones. To facilitate this, we will normalise each vector as we 
create it, and we will also maintain our list of vectors in order of decreasing 
weight. In fact, ( |34| ) tells us dim(V fc °), but we will build 25° as if we didn't know 
this. 



Recall that in §5.1, we created Z graded levels of B by repeated applications 
of the set of odd lowering generators to the set {|1}}. Here, each application of 
ESLGS generates an ungraded level, which we shall call Li, where Lq — 
We describe the process in the following algorithm: 

i := 

Lo :={|*i» 
while L, ^ 

L l+1 := (A(ESLGS) • L<) \ {0} 

normalise 

Li+i := Li + i \ L 

L := L U L i+1 

sort L by decreasing weight 
increment i 



Note that we have taken notational liberties by writing A(ESLGS) ■ Li. We 
do this as it is natural to apply functions to lists in Mathematica. 

Of particular interest here is that the evaluation of the algebra-module ac- 
tion in the tensor product case utilises the information encoded in the matrix 
elements of the underlying representation to determine when basis vectors are 
annihilated. That is, as different from we do not have to explicitly imple- 
ment the annihilation rules dictated by the KIMC. This observation carries over 
into the following subsections. 

13 We might just as well use ELGS, the full set of even lowering generators. The tradeoff 
is that whilst the coproducts of the nonsimple generators are more complex, there should be 
less levels to calculate. 
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6.4.1 Illustration: *B? for U q [gl(3\l)} 

We illustrate the process by constructing the 3 dimensional 25^ for i7 q [gZ(3|l)]. 
Recall that we determined: 

|*i> = ^o (? f |2®l)-g* |1®2», (76) 

hence L = {C^ (gf |2®1> -q% |1®2))}. We also have ESLGS = {E 2 U E 3 2 }, 
for which, from (E5a), we have the coproducts: 



A(E\) = {E 2 1 ®K*K1) + (K*K?®E 2 1 ) 
A(E 3 2 ) = (£ 3 2 ® K^Kl) + {k\k\ ® E 3 2 ). 



(77) 



Applying the operators of (77) to (|76|), we obtain two vectors: 
{0,eJ(3S|l®3)-g3|3<»l»}. 

We discard the 0, and find that the second vector is already normalised and 

i 

also not found in Lq, so we have L\ = {Cq (?p |1®3) — q% |3®1))}. As the 
components of Lq and L\ are independent, and indeed ordered, we have, at this 
stage: 

L = {Cl (tfi \2®l)-q^ |1®2}) ,cl (q9 |l®3}-g^ |3®1))}. 

Repeating this process on Li, we find L 2 = {cj (gf |4®1> -gf |1®4))}. As 

before, we normalise (again already OK) and install this vector in its rightful 
position in our collection L, checking first to see if we have already met it. We 
thus have, after ordering: 



L = 



Cl (q9 |2®1) |1®2)) 
Cl (g* |1®3) -q% |3®1)) 
Cl (q9 |4®1) -g^ |1®4)) 



Repeating again, we discover that L 3 = 0, so the process is completed. Our 
orthonormal basis 2$? is thus the above L. 

Note that, apart from some factoring, we have left this basis in the raw 
form that our Mathematica code yields. The human calculator, preferring 
symmetry, may wish to multiply some of the vectors by (the phase factor) — 1, 
but this is unnecessary for our purposes. In Appendix [b|, where our results are 
summarised, we have make some judicious changes of this nature for readability. 
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6.5 A nonorthogonal basis *Bfc for Vk 

At this stage, we have established the orthonormal basis 25° for V k °, and we 
wish to extend 25" to a basis 25 & for Vj,. This process involves two stages: 

• Firstly, we use the KIMC to construct a basis 25 & for Vk by the repeated 
actions of the odd lowering generators on the basis 25°, normalising and 
casting out repeated vectors as we go. This part of this process is detailed 



in this subsection: its appearance is similar to that of §6.4 



Unlike 25° however, the vectors of 25fc are not guaranteed to be orthogo- 
nal, as Vk in general contains some weight multiplicities. That is, distinct 
vectors of the same weights will generally appear, and these are usually 
nonorthogonal. To deal with this problem, we apply a Gram-Schmidt 
process to orthonormalise 25 fc into the final 25^. To optimise this, we 
preprocess 25 k by ordering its vectors by decreasing weight and then par- 
titioning it into weight equivalence classes. We then need only apply the 
Gram-Schmidt process to each equivalence class. The end result is 25fc, 
the desired orthonormal weight basis for Vk- This process is documented 



in §3.6 



Thus, we reproduce essentially the same algorithm as that used in §6.4, the 
essential differences being that the levels Li are now Z graded levels of Vk , and 
that we act with OLGS (the (full) set of odd lowering generators, see (f45|)), 
rather than with ESLGS^ 

i := 
Lo := 25° 
while Li ^ 

L l+1 := (A(OLGS) ■ L t ) \ {0} 

normalise Li + \ 

Li+i '■= Li + i \ L 

L := L U L l+l 

sort L by decreasing weight 
increment i 
®fe :=L 



Again, the elements of 25^ are only unique up to phase factors, and our 
code doesn't select these, so the final results contain various factors of —1 that 
a human calculator would quickly purge. Furthermore, the weight ordering 
covers some elegant symmetries of the generators. In the results presented in 
the Appendix, we make some judicious cosmetic changes for readability. 

14 In this case, if we try to only use OSLGS, the set of simple odd lowering generators, 
then we miss some of the vectors of each level, which are obtained by nonsimple odd lowering 
generators, i.e. products of simple odd lowering generators with even lowering generators. 
The combination of the use OSLGS and ELGS (or, the repeated use of ESLGS) to build 
each level would require more calculations. 
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6.5.1 Illustration: 03 1 for U q [gl(3\l)} 



We ill ustrate the results using <Bi for U q [gl(3\l)], continuing the example of 



.4.1. Here, 23i has 24 elements, sorted into equivalence classes of decreasing 
weight, and judiciously factored. To save space, we have written V = q + q, and 

v^v- 1 . 



_1 , a 

jg ( -7J~2~ |1®2> + q^I 3f{ 



C A 2,l ( A 1? <+7° +2 2®3) + ,= +2 |3C32>) + |1®5> - |5c31>)j 

V' f+9 5 |2tg3) + 7/3 |3®2> ] 



- |3«3> 



1 

c? 


^+,2+2 |3C85> - 


7J2+2 |5®3>j 






1 

7*"-' 

°0 


(-77 "3 1®4} + <jf 


|4gil)j 






1 


1/1 „, 1 

A 2 2 l [ A o C-5 + 2 


2®4> - , + 2 


1 

4®2)) + A 2 


(- 1®6> 4- 6®l))j 


1 

V 2 


[ -,2 |2®4> -7/2 


|4®2>^ 






1 

C? 


^ +s § + 3 |2C86> - 


7j3 + 3 |6®2)^ 






1 

°o 


i / 1 ,1 
A 2,l [ A < + ' 2 


3®4> + £J T 2 


1 

4®3>) + A.J 


(+ |1®7> - |T® 1>)J 


i 


/ 1 1 

+ g3 |3C34) + 7j3 


|4®3> j 







^+-A| ( + « 2 I 1 ® 8 ) - «^ 8g>l}")+A 2 ('+?2+ 1 |2®7> - ,2 +1 7«2>") 



'0 "1 "2,1 



V 2 Cj* 



+ a' (-7J 2 |3®6> 4- , 2 |6®3)J+a2 
+ 1 
+ 1 



'1 
1 

V 2 C 



1_ 1 

C l A 2 2 ,3 


/ 1 1 1 I \ 
Aj 2 (-f 2 5CS6) + q 2 |6®5>) + A J ( + + 1 |2(gS> + 7J" + 1 8® 2)) 






+ g § + 3 | 3[8 7> -ljf + 3 |7®3>^ 





°1 A 2 3 ( A l < + « 2 l 5 ® 7 > - ' 2 I 7 ® 5 )) + A 2 (-q" + 1 3®8) -7J"+1 8®3>) 
- |4®4> 



1 


'+3 + 3 


a , 1 

|4C86> - 7j3 T 3 


|6®4>^ 


1 

c 2 


-,§ + 3 


SL , 1 
4tg7) 4- 7f 2 T 2 


|7®4>j 
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6.6 Orthonormalising *Bfc into 23^ 

At this stage, the KIMC has yielded a nonorthogonal (although normalised!) q 
graded symmetric basis <8fc for Vk ■ With a view to constructing the projector 



Pfe onto Vk (see §7.1), we require dual bases for Vk and its dual V k *. From 



knowledge of 03 ^ there are two obvious ways to construct these dual bases: 

• We might continue to regard 03 & as our basis for Vk, and construct a 
(nonorthogonal) dual tensor product basis by the inversion of an overlap 
(i.e . metric) matrix. This process is described for <8i for [/ g [<7Z(3|l)] in 
fl3[| , although those authors don't actually implement it. 

• Alternatively, we can orthonormalise 23fc into Q3fc using a Gram-Schmidt 
process. The dual basis 23£ is then naturally orthonormal, and indeed 
trivial to write down. 

Implementation of both of these methods has shown that, apart from being 
substantially more efficient, the latter method yields more tractable and sym- 
metric results, so we choose it. Not only is it more elegant, but happily, the 
Gram-Schmidt process also allows us to maintain the q graded symmetry of the 
basis vectors. 

A basic principle in numerical computation is to never invert a matrix un- 
less absolutely necessary, as the process is both computationally inefficient and 
(numerically) unstable. The same comment about computational inefficiency 
certainly holds for the inversion of symbolic matrices. More seriously, for our 
current purposes, the inversion can bog down altogether due to difficulties in 
the simplification of algebraic expressions; a feature we might call 'symbolic 
instability'. 



6.6.1 Illustration: <Bi for U g [gl(3\l)] 



We illustrate the Gram-Schmidt process by continuing the example from §6.5.1 



that of 25i for U q [gl(3\l)]. Observe that the partitions of 23i include three of 
size 2 and one of size 3. To convert 23 1 to 23 1, we must orthogonalise each of 
those partitions. To illustrate, for the largest partition, we obtain the 3 vectors: 



°0 °1 2,1 



a| (+qi 1(88) - g§ 8®1)) + A J ( + qi + 1 2® 7) - ,§ + 1 7«2>) 
+ Aj ^-^f 3g>6) + ? f |6®3>) + A| + 1 |4®5) +7j^ + 1 |5®4>j 

A A 2 A 1 (~«^ I 1 ® 8 ) +9^ IS® 1>) + A S^ (-q% + 1 2®7) + + 1 |7»2>) 
1 / 3a , , 3a , , \ i / a a 



+ A2 



' 1^2,3 



A 2 A 2 A z 
A o J± ^ 



-„-2- 

(+q~2 |1®8> - <f~2 



f 1 |386> + qT + 1 |6®3> I + A2 I -q~Z 4 ® 5 > + q~S |5®4) ) 



31) + A 2 A-f 



[2®7> +q 2 



7X2} 



4 1 2 |3®6> + ,2 +1 6®3>) + A? (+7/2 |4®5> - ,! |5®4>) 



Observe that, again, our Mathematica code has done some nontrivial work 
in simplifying the algebraic expressions in q and a. This work would present a 
significant barrier for a human calculator. 
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7 Projectors and R matrices for V V 



7.1 Projectors Pk onto the 14 

At this stage, for each of the submodules Vk C V ® V, we have an orthonormal 
basis Q5fc = (l 1 ^?)}^™ ^ and, for each of their duals V k *, a corresponding dual 

basis W k = {(Vf \}fj^ iV "\ where <*f | = |*f) f . Using these dual bases, it is a 
simple matter to construct the projectors Pk : V <8> V — ► Vfc: 

dim(V fe ) 

p fc= E I I i 

J = l 

note that we must use (^2|) for the multiplication of tensor products. 

We now make a change of notation. As we did for the matrix elements in 



§5.2, we replace with the elementary matrix £ Mi^. We then use the 



notation ej* £ M^2mn to indicate the two dimensional matrix form of the usual 
elementary rank 4 tensor e* ® ef , obtained by inserting a copy of ef at each 
location of e'-. We find that the Pfe are in general quite sparse, that is, only a 
small fraction of their 2 components are nonzero. 

7.1.1 Illustration: P Q for U g [gl(3\l)] 

We illustrate using P for the case U q [gl(3\l)}, which has 125 (out of 2 12 = 4096) 
nonzero components. We present these components below, using horizontal lines 
to separate equivalence classes of symmetry. 

{ell} 



i (-{elleileltn ~ f^M\ a 2 [~4f 

{„23 „24 „341 r i f„27 „36 „45l 

e32,e42,e 43 l Ao_ I + |e 72 , e 63 , e 54 j 

{eM, e £}\ A eS}* 

+ {e?I, e 2 |,e? 4 7 }/' 9 \- {el?, e t?,e#f } 



f„32 42 431 J ! y " <| , f„23 p 24 341 
\ e 15i e 16> e 17/ J (,+ \ e 51i e 61i e 71/ , 

44 ^ 1+9^1, -e?l,+gef 8 /' 9 1 -<zeg, +e|?, -leg 

+g e l|,-efl,+ge!f/' 9 \+ge|I, -elf , + 9 c|f J 
-gel, +el, -gelfl _ a /— geg|, +eg, -gegl 
-?e| 6 7,+ef 7 ,- g efl/' 9 \~ge 7 |,+e^,-gelf/ 
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7.2 R matrices R m > n (u) and R m ' n 

We may now form the trigonometric R matrix R m,n (u) as a weighted sum of 
the projectors, where the weights are the eigenvalues of R m ' n (u) on the sub- 
modules. For the special case of our representations labeled A = (0 m |d„), these 
eigenvalues are actually known ]TT[ ] . The quantum R matrix is then the spectral 
limit R m ' n = liniM^oo R m - n (u). Again, what follows here strictly applies only 
to the to n case, but given the natural isomorphism between U q [gl(m\n)] and 
U q [gl{n\m)], this is unimportant. 

Again using the notation introduced in §|[ especially noting the definition of 
A 7 in (|64|), we have the following expression for R m ' n (u), normalised such that 
its 'first' component (i.e. the coefficient of e}}) is unity (for applications, other 
normalisations may be applicable, e.g. see 0): 

R m > n (U) = £ S 2A+A,, (U)P 2A+ A, , (78) 

[7] 

where, again, as in (|65[), the sum is over all allowable Young diagrams [7] and 
P 2A +A 7 is the projector onto the submodule V2A+A-. ■ Recalling that we intend 
7 = 7i. 72, ■ ■ • , 7r, the eigenvalue S 2 a+a 7 (u) is: 

. , -r-r -r-r \Oi + 7 — 7 + 77 L 

where, for the empty Young diagram case, we intend S2a(w) = 1. Note that we 
have substituted x = q~ 2u from the original x used in the multiplicative Yang- 
Baxter equations of p"l| ] : our Yang-Baxter equations are additive in variable u 
(seeQ. 

Then, as R m ' n is the spectral limit liuiu^^ R m,n (u), we have: 

Pr^ = j2^ 2A+Xy p 2A+Xy , (so) 

[7] 

where ^2A+a t (m) = liniu^oo S 2 a+a 7 (ti). Again, the coefficient of e{{ is unity. 
Evaluating this limit: 

ums 2A+A7 (u) = nn(-)? 2(a+r!) = n r^ 2 ^- 7 ^, 
j=i i=i j=i 

where we have applied the observation: lim-u^oo = —q +2X . Strictly, this 

requires that \q\ > 1, and this is perhaps not so sensible, as q is in some sense, a 
deformation parameter, so should be small. If, instead, we assume that \q\ < 1, 
then the limit becomes —q~ 2X , and the expressions for £ 2 a+a both above and 
below (in (pl|)), remain valid under the mapping q q. In the final analysis, 
this simply means that we obtain quantum R matrices related by q t— > q. As 
above, we are considering generic q (i.e. q is not a root of unity), and to ensure 
that, we can demand \q\ =/= 1. 
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These considerations aside, we obtain: 




intend £2 a = 1, in agreement with E,2a{u) = 1. 

Thus, we have two methods to compute R m ' n . Firstly, we may explicitly 
evaluate it as the spectral limit of R m ' n (u), itself computed by substituting ( f79|) 
into (|78|). Secondly, we may directly substitute (|8l]) into (|8(]), bypassing the 
construction of R m ' n (u) altogether. This method is of course less computa- 
tionally expensive, so in practice, we use it, but we have also implemented the 
former method, which is useful for checking consistency. 

7.2.1 Illustration: The R matrix decompositions for U q [gl(m\l)] 



For the special case n = 1, which we are interested in for physical applications, 
these simplify to: 



Again resurrecting the abusive notion of (|66|), we may replace the explicit ( f78| ) 
and @ with: 





k 



k 



rn 



m 





k=0 



where Sfe(tt) and £/. are the following eigenvalues and cf. (Tsj): 




here we of course again intend (u) 



1. 
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7.3 Yang— Baxter equations 

To be certain, R m,n (u) satisfies the following graded version of the (additive) 
trigonometric Yang-Baxter equation (TYBE): 

H IMM + M[c] +[ .]M + M^ tt) ^' R(u + v) ^ k{v) ^ 



where we have written [a] for [\a}]. The parity factors in (82) may be removed 
by the following transformation (e.g. see Q): 

K ab K U ) {-) K ab W> 

after which R(u) which satisfied ( p^ ) now satisfies the usual ungraded TYBE: 

x, ,c 6 ,c a 1 ~ , .b a 1 ^/ \b ' a ' ~, .c'a~,.cb 

i?( M ) b , c , i?(u + u) o , c i?(u) ab =R{v) a , b , R(u + v) ac , R(u) bc , (83) 

written in noncomponent form as: 

Ru(u)R 23 (u + v)R 12 (v) = R 23 (v)R 12 (u + v)R 23 (u). (84) 

In the spectral limit R = ]hn u — too R(u). t this of course becomes a quantum 
Yang-Baxter equation (QYBE): 

R12R23R12 = R 23 Rl 2 R 23 , (85) 

viz (R ® 7) (7 ®R)(R®I) = (7 <g (7? ® 7) (7 (g) , familiar as the braid relation 

OlO^CTi = CT2C1C2. 

Defining i?(it) = PR(u), where P is a permutation operator, yields a trigono- 
metric R matrix R(u) satisfying the following version of (p4|): 

i?i 2 (w)i?i 3 (u + v)R 23 {v) = R 23 (v)R 13 (u + v)R 12 {u). (86) 
This transformation amounts to the mapping: R{u) a ab = R(u) b ab . In compo- 



nent form, fl86| ) is more symmetric than (|83|): 
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7.4 An alternative construction of the quantum R matrix 

In and §0, we described the construction of trigonometric and quantum R ma- 
trices corresponding to the tensor products of representations of highest weight 
A, from explicit knowledge of the decomposition of the tensor product, and 
the eigenvalues of the R matrices on the subspaces of the decomposition. This 
method is limited to such situations where this data is known. 

An alternative approach to the construction of quantum R matrices sidesteps 
the construction of bases <8fc and projectors Pk altogether, instead using only 
the knowledge of the matrix elements. That is, say that we know the 'universal' 
(i.e. algebraic) form of the R matrix, i.e.: 



ft = ^Oi<g>&i, (87) 

i 

where a,i,bi <E U q [gl(m\n)\, for some hopefully finite sum over indices i. Then, 
for any particular representation of highest weight A, we may obtain a quan- 
tum R matrix tta®a(R) satisfying a parameter-free version of ( |S6| ) from 1Z, by 
simply replacing the Oj and bi with their matrix representations. In fact, for 
U q [gl(m\n)], in we find formulae for 1Z of the form (|S7|), so this method is 
feasible. 

Implementation of this approach has significant advantages over the current 
method; apart from being simpler, and greatly reducing computational effort, it 
is considerably more general in that it does not require knowledge of the tensor 
product decomposition. 

A substantial loss is that is does not yield explicit trigonometric R matrices, 
so we cannot use this method for physical applications - recall that our primary 
intended application is topological, being the construction of link invariants, for 
which we only require quantum R matrices. However, there is a further loss. An 
interesting alternative approach to constructing quantum R matrices involves 
constructing families of distinct but 'gauge equivalent' quantum R matrices, 
starting from a single trigonometric R matrix. (Details of investigations into 
this for the U q [gl(2\l)} case appear in [pp|.) Clearly, we also lose this alternative 
approach if we can't construct trigonometric R matrices. Furthermore, the 
method described in the present work incorporates the foreknowledge of the 
eigenvalues of the quantum and trigonometric R matrices. This knowledge is 
special to our representations, and it may be be used to assist analysis of the 



associated link invariants (again, see 1 20 ) . In the alternative construction, we 
have no such knowledge in general. 

The only barrier to implementing this method is that the results of Jl9| 
are presented in a somewhat abstract form which would require considerable 
modification before being directly useful in the framework described herein. 
This approach is outside the scope of this paper; it is left as a future project. 
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8 Implementation and results 



The entire process has been implemented as a suite of functions in the in- 
terpreted environment of Mathematica. The procedure to construct the R 
matrices requires several stages (viz the algorithms of §|| to §(?]) , and the several 
thousand lines of Mathematica code are broken down into functional units to 
achieve this. The code is available on request from the author. 



8.1 Data structure for the U q [gl(m\n)} generators 

A challenging issue in implementation is to find a consistent data structure to 
represent the algebra generators. This problem arises as the Cartan generators 
K a often appear exponentiated as , where N is not necessarily a positive 
integer. Whilst the unexponentiated K a (equivalently the E a a ) are a basis for 
the Cartan subalgebra, if we use them as our data structures, then we must 
deal with the problem of how to express and manipulate their exponentials. 
Thus it proves pragmatic to regard the more general 'generators' as the 
logical units for computation. Beyond this, our data structure must uniformly 
incorporate the non-Cartan E a b, which are never exponentiated. The follow- 
ing Mathematica pattern integrates these two disparate expressions into a 
coherent form: 

Generator [Uqgl [m, n] , a, b, N] , 

where m and n are the fixed m and n defining the algebra, and, if a ^ b (and N 
is fixed as 1), then we intend E a b, and if a = b, then we intend . 

Using this pattern, we are able to implement the collection of PBW commu- 



tators of §4.2 in only a few hundred lines of code. 



8.2 General comments 

Beyond the data structure, the following aspects of the code are specifically 
interesting as expositions of the use of Mathematica. 

1. Implementation of the PBW commutators to normal order strings of gen- 
erators uses the repeated application of Rules to find a fixed point. 

2. The code establishes and solves a system of linear equations to determine 
the parameters defining the highest weight vectors of the 14, where 
even the size of this system is not specified in advance (see §6.3). Doing 
this manually is particularly finicky due to the semantic complexity of the 
expressions involved. (In the literature, this process is specifically avoided 
whenever possible; e.g. [jl3| contains two different kludges.) 

3. Throughout, it has proven possible to maintain tensor product vectors 
and rank 4 tensor components in q graded symmetric combinations, which 
facilitates both the simplification of expressions and the presentation of the 
output. This has been achieved by applying rewrite rules to nonsymmetric 
expressions in a carefully controlled manner. 
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8.3 Limitations 



The computations are computationally inefficient, and this is due mostly to the 
fact that the algorithms used are direct, and not refined. Although there are no 
theoretical limits to m and n, computer storage and human patience mean that 
a current reasonable practical limit is ran ^ 4. Whilst the technically difficult 
translation of the interpreted Mathematica code into a compiled language 
would increase the speed of the computations enormously, storage requirements 
would still limit mn to perhaps 7 in the general case. 



8.4 Summary of results 

Both R" 1 ' 1 ^) and R" 1 ' 1 have been obtained for m — 1,2,3,4. Fixing m, these 
are rank 4 tensors, where the tensor indices run from 1 to 2 m (i.e. the dimen- 
sion of the underlying representation). Thus, they contain D m = 2 4m (albeit 
mostly zero) components. Where N' m and N m denote the number of nonzero 
components of R rn - 1 (u) and R™ 1 ' 1 respectively; as R m l is the spectral limit of 
R m ^(u), we find, as expected, that N m < N' m . We also find that N' m = 6 m 
(why?). Where s m = N m /D m denotes the sparsity of R" 1 ' 1 , we find that s m 
rapidly decreases with increasing m. Table presents this data for m — 1, 2, 3, 4. 



rn 


D m 


N' 




s m (%) 


1 


16 


6 


5 


31.3 


2 


256 


36 


26 


10.2 


3 


4096 


216 


139 


3.4 


4 


65536 


1296 


758 


1.1 



Table 4: Data for R matrix construction. 



For good measure, we also record the numbers of components of each of the 
?r; + 1 projectors for the above cases, finding them similar to N' m and N m . This 
data is included in Table |^, which also records submodule dimensions. 

m submodule dimensions projector sizes 

1 2,2 5,5 

2 4,8,4 25,34,25 

3 8,24,24,8 125,199,199,125 

4 16,64,96,64,16 625,1124,1254,1124,625 

Table 5: Dimensions of the submodules in the tensor product, and the number 
of nonzero components of the corresponding projectors. 
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Computer run times involved are listed in Table g. These show a rapid in- 
crease in cost with increasing m. This is accompanied by an exorbitant increase 
in storage required. 



in 





1 


2 


3 


4 


Representations 


0.10 


0.5 


2.8 


20 


TP submodule bases 


0.74 


7.4 


72.6 


1559 


TP projectors 


0.12 


4.0 


55.4 


877 




0.27 


6.5 


89.5 


1736 




0.04 


1.4 


18.3 


570 



Table 6: Run times for various calculations (in CPU seconds), using Mathe- 
MATICA 3 on a Sun ULTRA 60 computer. The timing for construction of the 
5 submodule bases for L/ 9 [gZ(4|l)] is the sum of the timings for the individual 
modules, viz: 1559 = 27 + 201 + 1061 + 218 + 53. 

Of these, the t/g[g>Z(l|l)] case can be done by hand in a few hours the 
complete J7 g [gZ(2|l)] case appears in my PhD thesis and took several weeks 
to do by hand; partial details (i.e. up to calculation of 3 out of 4 of the 58 for 
the [/g[<7Z(3|l)] case appear in |l|; whilst the £/ g [<7Z(4|l)] case is entirely new. 

Of interest is that R™ 1 ' 1 generally contains non-binomial irreducible polyno- 
mial factors as well as various q brackets. 

By direct substitution, we have been able to verify that each R m,1 (u) satisfies 
@0 and that each R m l satisfies @. 
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A Proofs of various lemmas 

A.l A commutation lemma — the proof of ( |T5| ) 

Lemma 2 

K a E c -q a E c K a , 

for any meaningful indices b, c, and any power N . 

Proof: 

Firstly, we show the following result: 

K a E b c = qf b - 6 ^E b c K a , (88) 

for any meaningful indices b ^ c, not just for the simple generators where we 
have \b — c| = 1. 

To see this, we first consider the case b > c, so that E b c is a lowering 
generator. We use induction on b — c to show the result, assuming ( pq ) for some 
b > c, where we know that it is true for c = 6 — 1 by (|l4|) . Then: 

K a E\^ K a (E b c E c c ^ 1 ~q c E c c „ 1 E b c ) 

= (K a E b c ) £ c c _i - g c (K a E c c ~i) E b c 

q f*- s °)E b c {K a E%^) - q c q { a"- 5 "- l] (E c c ^K a )E b c 

qi S Z- 5 ° ) E b Jf°- S ^ ) E%_ 1 K a - 

qJ a S *- 5 ^ ) EC c _ 1 qT^ ) E b c K a 

q W-r e +n-r e -i) ( EbcE c c ^ _ qcE c c _ lE b c ) K a 
ft) (st-s a c -i)p b K 

q a & c-l&a- 

Thus, (H|) is true for E b c -\ if it is true for E b c , and as it is true for c = b — 1, 
thus it is true for all b > c. The proof for raising generators follows by a trivial 
analogy. 

Next, observe that setting b = c in (p8|) also yields a true statement. That is, 



0) 



then states K a E b b — E b i } K a , which is equivalent to (12) when we replace 
K a with E a , and expand the exponential as a power series. 

Finally, replacement of q with q N in ( |88| ) and again using K a = q^ 1 ' B a 
allows us to write the more general statement in (|l5|). □ 

Having proven this, in retrospect, it appears that some use of the 'generalised 
Lusztig automorphisms' @, |2^| might facilitate a more elegant proof of the 
recursion in (|88|). 

A. 2 A coproduct lemma — the proof of (|28|) 
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Lemma 3 

A(E a b ) = E a b ®Kl S >Kl SS +K h a S »K^ S » ®E a b - 

S£ A <= (Kl St K^ St E a c ® E c b Kl Sl K\ S? ) , 

c 

for any valid indices a,b (including a — b), where the sum ranges over all c 
strictly between a and b, and is simply ignored if \a — b\ 1. 



Proof: 

Firstly, note that, for simple generators i.e. 
coproduct of (p7f): 



A{E% 



E a b ®Ki St Kl St +Kl Sl 



K 



1 

J h 



: — b\ ^ 1, the result is just the 



E a h 



as the sum is ignored. Specifically, it applies to the Cartan generators a = b, 
being A(E a a ) = E a a ®ld + ld® E a a , which is equivalent to (|25|c) once we make 
the identification K a = q(~^ 1 )e a . 

More generally, given (|27]), and the expansion of the nonsimple generators 
(||), the following straightforward induction on \a — b\ shows that the result 
follows for arbitrary nonsimple generators E a b . 

We first deal with the lowering case a > b, viz = +1. For the inductive 
step, we assume the truth of our statement for a — b = p, for some p > 1, and 
we show that this implies its truth for a — b = p + 1. Thus, beginning with: 



A(E a b 



-l) - 
(£1 



A (E a b E b b ^i — q b E b b -iE a b ) 

A (E a b) A (£Vi) - %A {E\_ 1 ) A (E a 



and expanding A (E a b ) using our inductive hypothesis and A [E b . 
definition (|25|a), we obtain: 



from the 



A(£Vi) = 

' +E a b E\- 1 ®KlKlKlKl_ 1 
+KZK*k}k}_ 1 ® E a b E\^ 



-KlK*E\^®E a b K*Kl 



-q b E\-iE a b ® K^Kl_ x K2Kl_ 1 
-q b K\iq_JilKl ®E\^E a b 
-q b E\_{KlKl 



K b 2 Kl 



rpa 



-E a b K b K, 



b-i 



Ka K b E b b _i 



-q b KlKl_ x E a b ® E\- Y KlKl 



a-l 

E 

2=b+l 



(Qc - q c ) ' 



+K 2 C K£ E a c E\^i ® E^KIKIK^KI^ 
-q b E\-iKlKlE a c ® kIkI^E^kIkI 

+KjKlE a c Kliq_ 1 ® E\K\~K\E h b-i 
-q b KlK^_^lK^E a c ® E h b ^E c b KlKl 



(h +t 2 +t 3 + U) 



a-l 

E 



A c (t 5 +t 6 ). 
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Note that in the above, no parity factors appear, as E a b and E b b _i cannot both 
simultaneously be odd. Next, we examine the terms 

(E a b E\^ - q b E b b ^E a b ) ® KIk\_ x ( l } £Vl <g> JQk\_ x 
K\Kl_ x ® (E a b E\. x - q b E b b ^E a b ) ( l ) K k a Kl_ x ® £Vi 

(kIkIe\^ ® e^kIku) (i - ash* ) = o 

{KlKt 1 E a b ®E b b ^K? l Kl) (qlql-q b ) 
-A b (kIk^E^ ® E h b _ x K2Kl) 

(k\k\e\E\- X ® E%KSkIkIkI_ 1 ) (l - g 6 gjg*) = 

k\k^ x E\ ® (E c b E\^ - q b E\^E c b ) k\k\ 
KlKl^E a c ® E c b ^K2K\. 
Thus, we have: 

A(£Vi) - E\_ x ®KlK\_ x +Kliq_ x ®E a b ^- 
A b (kIk^E^ ® E\^kIk^ - 

A c (K^K^E a c ® E\^KlKf) 

c=b+l 

= e\^ ® + # W_i ® ^ a 6 -i - 

a-l 

^ A c (KlK^_ x E a c ® E c b ^K%Kl) 

c=(6-l)+l 

Thus, the result is true for a — 6 = p + 1 if it is true for a — = p, hence it is 
true for all lowering generators £' a b, where a > b. 

Now observe that the definition of A for simple raising generators is ob- 
tained from that for simple lowering generators by the mapping q i— > q. Also, 
the definition of the nonsimple raising generators is obtained from that of the 
nonsimple lowering generators under the same mapping. Together, these defini- 
tions imply that the expression for A for the nonsimple raising generators may 
be obtained from that for A of the nonsimple lowering generators under that 
mapping q^q. 

□ 

It is also possible to deduce the coproducts of nonsimple generators via an 
entirely different approach, using the operator L(x) defined in |27fl , but we have 
not followed this up. 



h = 



5G 



B Explicit results for the £7jg/(3|l)] case 

Here, [i] = for i E {1; 5, 6, 7} and [i] = 1 for i 6 {2, 3, 4; 8}. Apart from a ll the 



notational conventions mentioned in the main text (recall A = q — g of §5.2.2 
and the Afj and Cij of j72] ) and ( f73| ) in §3.3), we add a couple more to condense 
the results. 

• We use the following notation as a shorthand for the q graded symmetric 
combination of tensor product vectors: 

q x ±\i®]) = q x \l®j)±q x \i®i) ■ 

This notation leads to such eyesores as "g° In these expressions, 

we shall choose i ^ j: observe that if j > i, we may replace q± \i®j) with 
±q ± x \j®i) = ±q x ± |i®i>. 

• To convert the graded R matrices into the equivalent ungraded objects, 
simply multiply all terms in boldface by —1. 

Having done this conversion the following notation is a convenient short- 
hand for the q graded symmetric combination of rank 4 tensors in the 
resultant ungraded R matrices: 

x ik A „x ik ki 

Q± e ji — 1 e ji ^ 9 e ij ■ 
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B.l Matrix elements of the U q [gl(3\l)] generators 

For completeness, here we present matrix elements 7i"(o,o,o I a)(X) f° r a ^ the 
[/ g LgZ(3|l)] generators X (including even the E a a , for comparison with the K a ). 

K{E\) = -ei-et-el-e% 
n(E 2 2 ) = -el-el-e 7 7 -el 
n(E 3 3 ) = -e\ -e 5 5 -el -e s 8 

7r(ifi) = e\ + e\ + e| + qe\ + e\ + qe% + qe 7 7 + qe% 
tt(K 2 ) = e\ + e 2 2 + qe\ + e\ + qe\ + e% + qe r 7 + qe% 
n(K 3 ) =e{ + qe\ + e\ + e\ + qe% + qe% + e 7 7 + qe% 

tt(E 4 4 ) = ae\ + (a + l){e\ + e\ + e\) + (a + 2)(ejj + e% + e?) + (a + 3)e| 
w(K 4 ) = q a e\+q a+ \4 + e§ + et) + q a+2 {4 + 4 + 4) + T + "4 
tt^s) = -e\ - e% 
n(E\) = -ei-4 
ir(E 2 3 ) = -e\ - ef 
tt(£ 3 2 ) = -el -e\ ^ 

tt(E\) = Ale\ + AM + A! 4 + A\e\ 

n{E\) = A\e\ + A\e\ + A\e\ + A|e? 
tt(E 1 3 ) = e h 7 -qe\ 
7T(E\) = 4-q4 

n(E 2 4 ) = qA*4 - A^ + qAM - A\e% 
7r(£ 4 2 ) = qA\e\- A\e\ + qA\e\ - A\e% ^ 
ir(E\) = -q 2 A\e\ - qAfe 2 - qA^e 3 7 + AJef 
n(E\) = -q 2 Al4 - qA\e% - qA\e\ + A*e* 

B.2 The basis <B for V = V( ,o,o|2a) 
The 8 vectors in this basis are: 

1 1(81) 

<7*g|{|l®2),|l®3),|l®4>} 

f4 9 + |i ® 5>+ 44 |2 ® 3> '] 

CcMf.i I A\q% |1®6) + Ah* |2®4> , > 
( A\q% 1 1® 7) + A\q\ |3®4) J 
cIclAl, [4 (<?| |3®6) - |2®7) - qf +1 |4®5>) - 4 5? |1®8>] 
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B.3 The basis «8i for Vy = V {0 ,o,-i 1 2a+i) 
The 24 vectors in this basis are: 

{|2®2),|3®3),|4®4)} 
c\ql {|1®2),|1®3),|1®4)} 
c\ ql + ^ {|2®5} , |3®5) , |2®6) , |3®7) , |4<g>7) , |4®6)} 
{Ajq'L\l®Z)+Ajt + + ^\2®2,)^ 

c!aI a I au°_ |i®6) + Ah+ + * I 2 ® 4 ) . ' 

[ Alq°_ |1®7) - |3®4> s 

f4g-|i®6>-4^ + ^|2®4}/ 

CMI.i < Ajg° |1®5) - A* |2®3> , ' 
{aU°_ \l®7)-A*q° + *\3®4) ( 

(4?? +1 |2®8>-4^|5®6),' 
CMI.s ^ Af^ +1 |3«>8> - A? ?! |5®7) , > 
( AJ^ +1 |4®8) - A\q\ j6®7) J 
^Ma,i [a| 9 J |1®8) + Ajgl +1 |2®7> - A ( f q f |3®6) - A,f + 1 |4®5>] 
CiAf Aj a [a| Algf |1®8) + A g f +1 |2®7) + A ig * +1 |3®6) + Ai«j f |4®5}] 
CiAf A| 3 [a* Ajq* |1®8> - A 2 g* +1 |2®7) - Aiq f +1 |3®6> + Ai<7* |4®5}] 

B.4 The basis <B 2 for \/ 2 = V (0 _i _i | 2q +2) 
The 24 vectors in this basis are: 

C{a\~a\ a [aIaI<iI +1 \1®8) + A qf +2 |2®7> + A iq f |3®6> - A l9 | +1 |4®5>] 
CiA? A| 3 [a| Ajg| +1 |1®8) - A 2 gf |2®7) - A^** 2 |3®6) - A 1? | +1 |4®5>] 
C?C|A| 3 [a| «jf +1 |1®8) - A|g| |2®7) + A|g| +1 |3®6> + A* g* +2 |4®5>] 

f A|^ +1 |1®5) - Af ql |2®3) , j 
£?4,i | A ( f^ +1 |1®6> - |2®4> , > 

[ 4^ +1 |1®7) - Af |3®4> J 

!A| g° j2®8) — Af <j" + ^ |5®6) ,1 
Af g° j3®8> - Ajt + +i |5®7> , > 
A|g° |4®8) - A? |6®7> J 
{A\q°_ |3®8}+a| 9 + + J |5®7),1 
4,3 \ A\ q°_ |4®8) + Af <£ + 5 |6®7> , > 
( A?g° |2®8) + Af g" +i |5®6> J 
+ 5 {|2®6} , |2®5) , |3®5) , |3®7) , |4®6) , |4®7}} 

clqf +1 |j5®8) , |7®8) , |6®8)} 
{|7®7},|6®6>,|5®5}} 
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B.5 The basis <B 3 for V 3 = V(_i _i _i 1 2a +3) 
The 8 vectors in this basis are: 

c\c\~a\ z [Alqf +A 1 1® 8) - A\ (q f+2 |2<87> - g f +1 |3<8>6) + g f |4®5) 

A\ q1+ 2 |2(8>8) + Af |5(86) , 

^1^1,3 { Afq1+ 2 |3®8> + |5®7) , 
Afq%+ 2 |4®8) + Af |6®7) 

' i 2L + 1 



B.6 The trigonometric R matrix R 3,1 (u) 

For the listing of the components of ^ 3,1 (u), we invoke a little more notation: 
Sf = [a + i ± 

f/f = [u-i] q z , where z£ {1,2}, 
and i E {0, 1, 2}. With this, R 3,1 (u) has 216 nonzero components: 



Q+ Q+ Q+ Q+ Q+ Q+ 

1 /„H\ /„22 33 „441 "->() °1 r 55 66 771 °0 °1 °2 r„881 

l{en), — {e 22 ,e 33 ,e 44 \ , 7^7^ {e 55 , c 66 , e 77 j , c _ c _ c _ {e 88 ) 



S 5*! 



^0 ^1 ^2 



A, J<f {eg, e}|, e\i} \ A 2 S+ 5+ /?" {eg, e|l, el} 

Sjf 1 9" {elLelLet}} J ' s^fS^ \? u {eg, eg, eg} 

_i / /1 («) { ei , el! , el\ } \ 1 / h («) { e?g , e?g , egg } 

A^Sf {ei,ei,ei} J ' A'SoS^S; \h(q) {eg?,eg?,egg} 

AoAx / g 2 "{ei;,e^,e 7 i}\ 4^ fg 2 " {e||, e||, eg} 1 

^ ' {e&e&ei?}/' S^fsf \g 2 " {eg, el, el} J 
f/3(9){ei?},/4(g){ei},/ s ( g ){etl} 
/3©{e 7 |},/4®{e|I},/5@{el|} 



s -5r u 2u 



A 2 S S 2 

22_o J 1 \ e 52, 653, e 62 , e 64 , e 73 , e 74 | 



,26 „46 „3 7 „471 

°— e^ 7 , e 47 } 



A A 1 A 2 



Hi} 



U f+l{e^ 2 ,eM,eIt} 

, [ 21 31 411 

-1 {ei2,ei 3 ,ei 4 } 



S S 1 



s s x 

U0U1 
&o "1 ^2 



,,15 16 17 
e 51i e 61i e 71 
„51 „61 71 
e 15> e 16: e 17 

„23 „24 „34 
e 32) e 42i e 43 
„32 „42 „43 
e 23) e 24) e 34 



U0S+S+ f+l{ei 8 7 ,ej,e|l} 
SoS-S 2 \-l{e?Lel,el|} 

U0U1S+ 

"0 ^1 ^2 
^S + 



^0 ^1 ^2 
+l{el|,ei,e 2 72 7 }\ U £% 

1 r„45 „36 „2 7 1 

-1 }e 5 4,e 6 3,e 7 2} 



„84 83 82 
e 48, e 38> e 28 
„48 „38 „28 
^84, e83) e 82 

„76 „75 „65 
e 67) 657, e 56 
„67 57 56 
e 76) e 75i e 65 



U0U1U2 

S S 1 S 2 



+l{e 8 !} 
-l{e!I} 



1 / p 52 53 62 64 73 741 - 
1 \ e 25i e 35) e 26i e 46) c 37) e 47/ 
„_ „_ "1 | -i r 25 „35 „26 „46 „37 „471 

s o S i l +1 { e 52,c 5 3,e 6 2,e 6 4,e 7 3,e 7 4} 
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A§ Al Up 
S S 1 



1 ! „51 „61 71 1 1 \ 

-1 {e 3 2,e42,e4 3 j I 
+l{e|?,e|?,e|fr 

. -, f„51 61 71 1 

-l{e!?,eM} 

-l{e?|,e?tef 7 j 
-lfe&eJS.eJS 
+l{ e ?§, e «ef 7 



Sg 



+1 {eg, 
-l{ei, 
-Ijell, 
+1 {ell, 
-Ijell, 
+l{el, 
+l{ei, 



-l{r 28 



„75 „76l 

e 75, e 76/ 
57 671 
e 83, e 84 / 
„83 841 
e 57, e 67/ 

e%,e%} 

p 38 481 

e 3 7 l,ell} 
e 3 7 l,e|I} 



I ■• I • I •/< .J.v< | I -p /+! { e 36 5 ' 




S S 1 S 2 



where: 



Mq) 
Mq) 
Mq) 

Mq) 

Mq) 
Mq) 



-2q+(q 1+2a + q 1+2a )-q 2u (q-q) 
-2q+(f +2a +q i+2a )+q 2u (q-q) 



-tm 2 - (q 2+2a + q 2+2a ) + q 2u (q 2 ~tf)) 

g u (-2 + (q 2+2a + q 2+2a ) - (q 2 - 2u + q 2 - 2u ) + (q 2u + q 2u )) 
—ut n 2 . / 2+2a , —2+2a\ , 2u / 2 —2\\ 

q (-2g +{q +q ) + q (q - q )) 

2 + 2a 



q(q + q)-(q a + q 



.2+2^ 

.2+2a\ 



2u/ 2 

7 (q 

2u-l 



(«-«)■ 
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B.7 The quantum R matrix R 3 ' 1 
R 3 ' 1 has 139 nonzero components: 

1 iV 1 ! n 2a So 22 f> 33 ^ 44 \ r, 4a + 2 i^ 55 ^ 66 f> 77 l ^ 6a + 6 / < = 88 \ 

J-i e ii/, — ? \ e 22, e 33 , e 44 , j- , g te 55 , e 66 , e 77 ) , —q \ e ssj 

A„ a A So 21 o 31 c 41 \ A„ 5ct+4 4 Je 87 o 86 c 85 \ 

-Ag A {e 2 i,e 3 i,e 4 i) , -Ag A 2 (e 87 , e 86 , e 85 j- 

A 2a+l / 32 42 431 A„ 4a + 3 /„ 7 <> „75 „651 

Ag {e 32 ,e 42 ,e 4 3j , -Ag |e 76 , e 75 , e 65 j 
A 2 g 2 " +1 A,.4i {eg.egi.eB} , —A 2 q 4a+3 A 1 A 2 {elt,ei,ei} 
-A 2 g 3a+2 Ai {ei} , -Ag 3Q+3 Ai(g 2 - g 2 ) {e 7 l} 

A 3t>+l . r 52 53 „62 64 73 741 A 3 „ 3q + 3 4 4 A ! „ 8 1 1 

Ag Ai {e 52 ,e 53 ,e 62 ,e 64 ,e 73 ,e 74 j , -A g A)Ai,4 2 |e 8 i| 

-. r 12 ,13 141-1 |- 1 f„78 „68 58 

-1 {e 2 i,e 3 i,e 41 j I 5a+4 I -1 |e 87 , e 86 , e 85 

+l{e?2, e ?ieli}/' q l+l{e 8 78 7 ,ei,el| 

(e 15 P 16 f> 17 ~l fV 84 P 83 <= 82 

2a I e 51> e 61> e 71 I 4 "+2 |e 48 ,e 38 ,e 28 



9 S „B1 61 71 f , -5 S 48 38 28 

\ei5i e 16> e 17 J (, e 84i ^83) e 82 

C-23 -24 34-l ( 76 .75 .65 

2a+i I e 32 , c 42 , e 43 I 4a+3 Je 67 ,e 57 ,e 56 

-9 S „32 42 43 f , 1 S 67 57 56 

^ e 23) c 24) e 34 J 1^76) e 75i e 65 

,62 53 73 „64 741-1 |- 1 / „72 63 



3a+l 



1 ;„52 62 53 „73 „64 „741 -i ( -, / „72 63 „541~i ( i f„181 

1 { e 25 ) e 26 i e 35 > e 37 > e 46 j c 47 J I 3a+2 I ~1 { e 27 , e 36 , e 45 j I ■ ia \-l{e 81 \ 



1 ill f„25 „26 „35 „37 46 „471 ( , <1 i . r 27 36 „45 \ f il i . r 81 

I +1 {e 52 ,e 62 ,e 53 ,e 73 ,e 64 ,e 74 } I I +1 |e 72 , c 63 , e 54 } I l+l{e 18 



Ag A Ai|g 2 {_ l{ei 3 , ei 4 ;e 34j|,^ ( +l{e 32 e 42 e 43 } 

A 4 a+ 3, A Li f+l{ei,eILelI}l i f -1 {eg, eg, ell} \\ 

Ag Al^ 2 |g 2 | _ x {e|2; e|3j ,84} } , +1 {ei) £ 84 } ) ) 

^ +i {-{-;ii:J!}-« 4 {"il!}} 
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